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ABSTRACT 


RADAR  BASED  CASE  STUDY  OF  A  NORTHEAST  COLORADO  SNOWSTORM 

Radar  and  convoitional  meteorological  data  collected  during  the  19-20  January 
1991  winter  storm  over  Northeast  Colorado  as  part  of  WISP-91  (Winter  Icing  and  Storm 
Project)  were  analyzed  to  describe  observations  of  snow  bands  and  possible  band 
formation  mechanisms.  A  review  of  the  synoptic  situation  led  to  classification  of  this  sttxm 
as  an  anticyclonic  type  of  upslope  storm.  Single  and  dual-Doppler  radar  analyses  showed 
well  organized  mesoscale  snow  bands  oriented  west-southwest  to  east-northeast,  paraUel  to 
the  800-4(X)  mb  thermal  wind,  formed  during  upslope  conditions  several  hours  after 
passage  of  an  arctic  cold  front  over  the  research  area.  Embedded  within  the  bands  were 
precipitation  cores  that  propagated  along  the  band  with  the  upper  level  winds  and  nearly 
statKHiaty  (band  relative)  trapped  mountain  waves.  The  evolution  and  structure  of 
snowband  kinematic  pr(q)erties  were  studied  by  analysis  of  horizontal  and  vertical  velocity 
profiles  and  reflectivity  distributicms  using  PPI  scans  and  Extended  Velocity  Azimuth 
Display  (EVAD)  techniques.  Various  mechanisms  were  examined  for  their  potential  role  in 
formation  and  maintenance  of  snow  bands.  These  mechanisms  included:  conditional 
instability.  Conditional  Symmetric  Instability  (C!SI),  ducted  and  internal  gravity  waves,  jet 
streak  circulations,  seeder-feeder  cloud  structures,  and  terrain  induced  phentnnena  such  as 
tnqrped  mountain  waves,  low  level  barrier  jets,  and  cold  air  damming. 
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Chapter  1 
INTRODUCTION 


Northeast  Colorado  winter  storms  bring  a  wide  range  of  snowfall  amounts  with 
accompanying  beneficial  to  detrimental  ectmomic  and  social  impacts.  The  effects  of  a 
heavy  snowfall  are  far  reaching  and  may  have  national  in  addition  to  local  impacts  when 
the  role  of  Denver  -  Stapleton  Airport  (and  the  new  Denver  Airport  when  completed)  as 
a  major  airline  hub  in  the  nation’s  air  transportadon  system  is  considered  (Schlatter  et  al. 
1983,  Reinldng  and  Boatman  1986,  Dunn  1988).  Though  the  consequences  of  a  heavy 
snowfall  may  be  far  reaching,  the  actual  distribution  of  snow  can  be  highly  localized. 
For  exanple,  during  the  C!hristmas  Blizzard  of  1982  snow  depths  varied  between  over 
76  cm  (30  in)  in  the  Etenver  Metropolitan  area  to  less  than  5  cm  (2  in)  over  a  distance  of 
less  titan  60  km  to  the  north  (Schlatter  et  al.  1983).  The  orientation  of  low  level  winds 
with  respect  to  topographical  features  as  well  as  precipitation  bands  are  significant 
factors  in  snowfall  distribution  and  amounts  (Wesley  1991,  Wesley  and  Pielke  1990). 

In  view  of  the  important  nature  of  winter  storms  in  the  region,  the  Winter  Icing 
and  Storms  Project  (WISP)  was  initiated  in  1990.  In  this  thesis,  data  from  the 
WISP-91  (15  January  through  31  March,  1991)  experiment  were  used  to  examine  a 
snowstorm  that  occurred  19-20  January  1991  with  the  goals  of  discussing  observations 
of  siK>w  band  activity  and  evaluate  fonnation  and  maintenance  mechanisms.  Chapter  2 
provides  the  synoptic  and  mesoscale  meteorological  overview  incorporating  the  many 
assets  deployed  during  WISP-91 .  Radar  data  and  their  analyses,  important  elements  in 
obsovation  of  storm  structure  and  snow  bands,  is  the  subject  of  Chapter  3.  Possible 
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mechanisms  for  the  creation  and  maintenance  of  snow  bands  and  a  conceptual  model  are 
discussed  in  Chapt^  4.  Chapter  5  presents  conclusions  and  recommendations  for 
future  research. 

1.1  General  Background 

Until  recently  intense  field  experiments  of  winter  storms  in  the  United  States 
have  been  concentrated  over  the  coastal  areas  of  the  Pacific  Northwest,  and  Eastern 
Seaboard  with  only  a  small  number  of  experiments,  limited  in  scope,  conducted  over  the 
mid-ccmtinent  (Ramamurthy  et  al.  1991).  Hobbs  et  al.  (1980)  described  the  objectives, 
facilities,  and  q)erations  of  the  CYCLonic  Extratropical  Storms  Project  (CYCLES), 
which  took  place  over  coastal  Washington  State  and  adjacent  ocean  areas.  Numerous 
ptq)ers  rni  a  wide  variety  of  topics  related  to  winter  storms,  including  storm  kinematics 
and  miciophysics,  were  published  under  the  umbrella  of  CYCLES  research.  Similarly, 
Dirks  et  al.  (1988)  and  Raman  and  Riordan  (1988)  discussed  the  Genesis  of  Atlantic 
Lows  Experiment  (GALE)  occurring  off  the  Mid-Atlantic  States. 

literature  on  major  mid-continent  winter  storm  research  projects  is  more  limited. 
Ramamurthy  et  al.  (1991)  described  goals  and  objectives  of  the  University  of  Illinois 
Winter  Precipitation  Program  (UNIWIPP)  and  listed  other  important  studies  on  mid- 
ccmtinent  winter  storm  research  including  works  by  C^arbone  and  Bohne  (1975),  Bdine 
(1975),  Heymesfield  (1979),  Byrd  (1989),  Moore  and  Blakeley  (1988),  and  Agee  and 
Gilbert  (1989).  These  authors  discussed  winter  storm  activity,  including  kinematics  and 
mictophysics,  in  the  Great  Lakes  Region  and  relatively  flat  land  areas  of  the  central 
plains  of  the  United  States. 

Winter  field  projects  conducted  over  the  inteimountain  western  United  States 
have  been  even  more  limited  in  scope,  and  in  general  were  concerned  with  gaining  an 
understanding  of  the  processes  that  would  lead  to  precipitation  enhancement  through 
weadier  modification  (Rasmussen  et  al.,  1992).  The  Colorado  River  Basin  Pilot  Project 
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(Marwitz  1980),  conducted  during  the  five  winter  seasons  between  1970-75,  was 
designed  to  evaluate  precipitation  enhancement  potential  over  the  San  Juan  Mountains  of 
southern  Colorado.  The  Sierra  Cooperative  Pilot  Project  (Reynolds  and  Dennis  1986) 
was  a  weather  modification  experiment  performed  fiom  1976  to  the  late  1980’s  over  the 
American  River  Basin  in  the  Sierra  Nevada  of  California.  The  Colcxado  Orographic 
Seeding  Experiment  (Raub^  et  al.  1986)  studied  winter  clouds  over  the  Park  Range  of 
nothem  Coloado  with  the  goal  of  developing  a  weather  modification  hypothesis.  The 
State  of  Utah  Division  of  Water  Resources  and  the  National  Oceanic  and  Atmospheric 
Administration  (NOAA)  Tushar  Mountains  Cooperative  Weather  Modification  program 
(Long  et  al.  1990,  Sassen  et  al.  1990)  is  an  ongoing  study  of  weather  systems  over 
mountainous  regions  of  southern  Utah  designed  to  gain  more  insight  into  conditions 
favorable  for  weather  modification.  All  of  the  above  studies  were  focused  on  research 
into  conditions  and  weather  systems  favorable  for  precipitation  enhancement  Since  the 
storms  most  likely  to  fall  into  this  category  would  have  Pacific  origins,  these  projects 
were  located  on  the  western  side  of  the  ClcHitinental  Divide  thus  did  not  address 
conditions  over  eastern  Colorado. 

De^ite  the  lack  of  large  scale  field  projects  over  eastern  Colorado  prior  to 
WISP-91,  research  has  been  ctmducted  by  several  investigators  on  winter  storms  over 
northeast  Colorado.  Wesley  (1991)  provides  a  good  summary  of  significant  research 
(Ml  winter  sttMms  over  the  region.  Briefly,  winter  storms  generally  fall  into  two 
categories,  deep  cyclonic  and  shallow  anticyclonic  (Reinking  and  Boatman,  1986). 

Deq>  cycltMiic  systems  encompass  easterly  flow  to  levels  above  700  mb  in  response  to  a 
deep,  stacked  low  pressure  area  typically  centered  over  the  four  comers  region  of  the 
southwest  United  States.  Anticyclonic  systems  are  characterized  by  shallow  easteriy 
upslope  flow  over  the  eastern  plains  of  Colorado  (x;cuning  in  the  anticyclonic  flow 
asscxtiated  with  a  post-fiontal  arctic  high  pressure  area.  Hybiid  systems  encompassing 
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features  of  both  may  also  occur  (i.e.  andcyclcmic  easterly  surface  flow,  low  pressure 
and  cyclonic  easterly  flow  aloft;  Wesley  1991).  Alternatively,  climatological  studies  of 
heavy  snow  over  eastern  CoIch^o  in  progress  at  Colorado  State  University  (CSU) 
suggest  this  classification  scheme  might  need  modification.  The  location  of  a  surface 
low  pressure  area  in  the  vicinity  of  the  Colorado-New  Mexico  border,  which  is  a 
common  occurrence  in  all  cases,  is  a  significant  factor  in  snowfall  distribution  and  is  a 
better  descriptor  of  weather  type,  ratha*  than  the  cyclonic  and  anticyclcmic  classification 
described  above  (personal  communication.  Prof.  T.  McKee).  For  the  purposes  of  this 
paper  however,  the  classification  scheme  described  by  Reinldng  and  Boatman  (1986) 
will  be  used. 

1.2  WISP-91  and  Data  Sources 

The  Winter  Icing  and  Storms  Proj«:t  in  1991(WISP-91)  was  the  first  intensive, 
highly  instrumented  field  experiment  to  investigate  winter  storms  and  aircraft  icing  over 
northeast  Colorado.  The  project  evolved  frcnn  a  need  to  conduct  aircraft  icing  research 
and  a  genoal  desire  for  winter  storm  study  the  research  community  to  an  experiment 

with  these  two  broad  goals:  "1)  to  improve  our  understanding  of  the  processes  involved 
in  die  production  and  depletion  of  supercooled  liquid  water  content  in  winter  storms  and 
2)  to  improve  forecasts  ofaircrcft  icing  in  winter  storms"  (Rasmussen  and  Politovich 
1990).  Rasmussen  et  al.  (1992)  provide  an  overview  of  the  experiment  goals, 
operations,  organization,  participation,  and  instrument  platforms.  Figure  1.1  shows  die 
location  of  facilities  including  the  Colorado  State  University  (CSU)  CHILL  and 
National  (Dceanic  and  Atmospheric  Administration  (NOAA)  Mile  High  (MHR)  radars, 
forming  a  dual-Doppler  pair  used  during  the  experiment,  and  for  this  research. 
Appendix  A  lists  station  location  infOTmation  including  station  callsign,  name,  latitude. 
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longitude,  and  elevation.  More  details  on  data  collection  platforms  are  available  in  the 
WISP-91  Data  Catalog  (NCAR  1991). 

Surface  data  were  collected  every  minute  at  the  39  National  Center  for 
Atmospheric  Research  (NCAR)  Portable  Automated  Mesonet  (PAM)  stations  and  every 
5  minutes  at  the  NOAA  Prototype  Regional  Operational  Forecast  System  (PROFS) 
surface  mesonet  staticms  in  addition  to  the  hourly  and  special  observations  taken  at 
NaticMud  Weather  Service  (NWS),  Department  of  Defense,  and  Federal  Aviation 
Administraticm  weather  stations.  Complementing  the  twice  daily  NWS  rawinsonde 
network  observations  were  those  taken  at  3  hour  intervals  at  the  NCAR  Cross-chain 
Loran  Atmospheric  Sounding  System  (CLASS)  sites  between  1800  UTC  19  January 
and  21(X)  UTC  20  January.  Radar  volume  scans  were  completed  about  every  10 
minutes.  V^nd  profiler  data  were  available  every  20  minutes.  Given  the  large  variety  of 
meteorological  sensors  used  during  the  experiment  and  the  volume  of  data  obtained  for 
this  case  study,  many  analyses  and  display  t^hniques  were  used  to  evaluate  the  data. 
Appendix  B  describes  some  of  those  techniques  and  software  used  or  developed. 

Participation  in  the  data  collection  phase  of  WISP-91  at  the  recently  installed 
CSU  CHILL  10  cm  Doppler  radar  just  north  of  the  Greeley  airport,  offered  the  chance 
to  view  each  case  and  make  an  immediate  estimation  of  its  suitability  for  further  study. 
The  case  of  19-20  January  1991  was  selected  on  the  basis  of  data  quality  and  the  highly 
banded  structure  of  the  observed  echoes.  Using  the  Boatman  and  Reinking  (1986) 
scheme,  this  case  is  best  classified  an  anticyclonic  upslope  event  since  precipitatitxi  did 
not  begin  until  well  after  passage  of  the  surface  cold  front  where  the  region  was  under 
the  influence  of  an  arctic  high  pressure  area.  Chapter  2  provides  a  detailed 
meteorological  overview  of  the  storm  using  the  WISP-91  data  set  collected  for  those 
days. 

1.3  Northeast  Colorado  Topography  and  Climate  Controls 


5 


The  WISP  91  data  set  provided  a  unique  opportunity  to  examine  in  detail  the  role 
of  terrain  features  on  winter  storms  over  northeast  Colorado.  Previous  studies  by  Dunn 
(1988),  Wesley  and  Pielke  (1990),  Wesley  (1991),  Wesley  et  al.  (1990),  and  Toth 
(1987)  made  use  of  the  NOAA  Prototype  Regional  Operational  Forecast  System 
(PROFS)  surface  mesonet  which  has  greater  resolution  than  the  standard  airport 
observing  networks.  As  Figure  1.1  shows,  the  addition  of  39  PAM  stations  extended 
the  netw(»k  and  filled  gaps  in  the  coverage  with  respect  to  the  complex  terrain  found  in 
the  region. 

The  rt^id  terrain  rise  of  the  Rocky  Mountain  Front  Range  coupled  with  arctic 
ainnass  outbreaks  over  the  central  United  States  create  a  climatology  unique  from  that  of 
the  rest  of  the  United  States.  Figure  1.2  is  a  topographical  map  of  eastern  Colorado 
showing  the  con^licated  terrain  variations  frcnn  the  eastern  plains  to  the  Continental 
Divide.  Several  features  merit  notice:  the  majOT  west-east  extensions  of  higher  terrain 
embodied  by  the  Cheyenne  Ridge  near  the  Wyoming  border  and  the  Palmer  Divide 
south  of  Denver,  and  the  rapid  increase  of  terrain  elevations  found  along  the  Fixxit 
Range.  The  two  west-east  divides  separate  three  major  river  drainage  basins.  The 
Nortii  Platte  River  flows  north  of  the  Cheyenne  ridge,  the  South  Platte  River  Basin  lies 
betweoi  the  Cheyenne  Ridge  and  Palmer  Divide,  while  the  Arkansas  River  runs  south 
of  tile  Palmer  Divide.  These  extensions  of  elevated  terrain  significantly  modify 
predpitatitm  distribution  over  eastern  Colorado.  The  Christmas  Blizzard  of  1982  is  an 
excellent  exanqile  of  how  these  extensions  can  affect  snowfall  over  the  region  (Schlatter 
etal.  1983). 

Qimate  controls  for  tiie  eastern  plains  and  Front  Range  of  Oilotado  differ 
rignificantiy  from  those  governing  the  Piedmont  on  the  eastern  approach  to  tiie 
Appalachum  Mountains  (Lilly  1981).  Foremost  in  these  differences  is  the  lack  of  a 
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relatively  warm  ocean  and  maritime  airmass  to  the  east  Continental  arctic  air  can  plunge 
southward  out  of  Canada  into  the  central  United  States  unhampered  by  terrain  barriers. 
Colorado  terrain  is  significantly  higher  both  over  the  plains  to  the  east  and  the  mountains 
to  the  west  dian  ccxnpaiable  regions  of  the  Appalachians  and  adjacent  Piedmont.  In 
general,  snow  over  eastern  Colorado  occurs  after  a  cold  front  passes  or  development  of 
low  pressure  near  the  Colorado-New  Mexico  border  and  requires  easterly  flow  to  lift 
nooist  air  having  origins  in  the  Gulf  of  Mexico  against  the  Ront  Range  of  the  Rocky 
Mountains.  The  depdi  of  the  easterlies  combined  with  the  location  and  speed  of  storm 
scale  dynamics,  with  respect  to  terrain,  governs  the  location,  duration,  and  amount  of 
snow  (Reinking  and  Boatman  1986).  However,  there  are  cases  where  strong,  moist 
westerly  flow  aloft  over  a  deep,  extremely  cold  arctic  airmass  with  weak  easterly  surface 
flow  can  cause  prolonged  heavy  snow,  contrary  to  the  popularly  held  opinion  that  heavy 
snows  result  frimi  deep  easterly  flow  (Wesley  et  al.  1990). 

1.4  Precipitation  Bands 

Precipitation  bands  are  a  commonly  observed  mesoscale  feature  of  nearly  all 
winter  storms,  including  those  of  northeast  Colorado  (Rasmussen  et  al.  1992). 
Extensive  studies  of  precipitation  bands  were  conducted  during  (CYCLES.  Six  types  of 
lainbands  were  identified  and  classified  as:  warm  frontal,  warm  sector,  wide  cold 
frontal,  narrow  cold  frontal,  prefrontal  cold  surge,  and  post  frontal  (Matejka  et  al. 

1980).  Ptefiontal  cold  surge  bands  need  not  be  considered  here  since  snow  bands 
discussed  later  occurred  after  the  cold  front  passage.  Studies  of  precipitation  bands 
woe  also  published  as  part  of  GALE  (Emanuel  1988)  and  the  UNIWIPP  experiment 
(ShieUs  et  al.  1991).  Discussion  of  precipitaticxi  bands  over  northeast  Colorado  has 
been  limited  to  a  few  studies  described  below. 


7 


Hobbs  et  al.  (1980)  discussed  the  organization  of  cold  front  precipitation  bands. 
They  identified  three  types  of  bands,  the  warm  sector  (or  prefrontal)  band,  wide  cold- 
frontal  bands,  tens  of  kilometers  in  width  straddling  the  front,  and  narrow  cold-frontal 
bands,  ^Dproximately  4  km  wide,  located  at  the  surface  frontal  boundary.  As  before, 
the  warm  sector  band  is  of  no  interest  here  since  this  study  focused  on  the  period  after 
cold  front  passage.  AU  wide  cold-frontal  bands  observed  by  Hobbs  et  al.  (1980)  were 
oriented  parallel  to  the  front  and  moved  in  the  same  direction  as  the  front  though  at  a 
greater  speed.  Occasionally  wide  cold-frontal  bands  were  observed  to  overtake  the 
surface  front 

The  mesoscale  and  microscale  organization  and  structure  of  narrow  cold-frontal 
rain  bands  observed  during  CYCLES  were  described  in  detail  by  Hobbs  and  Persson 
(1982).  They  found  the  band  was  comprised  of  heavy  precipitation  cores,  ellipsoidal  in 
shqie,  with  an  average  orientation  angle  of  29*  between  the  synoptic  scale  cold  front 
and  the  long  axis  of  the  precipitation  core,  and  gap  regions  of  lighter  precipitation.  The 
cores  moved  with  the  speed  of  the  front  but  also  had  a  component  of  velocity  parallel  to 
the  frrxital  boundary.  Leading  edges  of  core  regions  were  also  marked  by  strong  low 
level  convergence  and  upward  vertical  motions  on  the  order  of  1- 10  ms'^  Modeling 
studies  by  Rutledge  and  Hobbs  (1984)  indicate  that  rapid  growth  of  ice  crystals  by 
liming  in  the  strong  updraft  of  the  narrow  cold-frontal  rainband  is  responsible  for 
precipitation  (which  can  exceed  100  mm  h*^). 

Wave-like  warm-frontal  rain  bands  observed  during  CYCUES  were  found  in 
advance  of  and  parallel  to  the  warm  front  (Herzegh  and  Hobbs  1980).  Typically,  these 
bands  are  die  result  of  a  seeder-feeder  mechanism  where  ice  crystals  fall  from  shallow, 
upper  level  convective  cells  located  in  potentially  unstable  air  (dOw/dz  <  0,  where  9w  is 
wet-bulb  potential  temperature,  and  z  is  height),  and  scavenge  cloud  water  in  lower 
level  stable  cloud  layers  (Matejka  et  al.  1980,  Houze  et  al.  1981).  Simulations  by 
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Rutledge  and  HoW)s  (1983)  of  two  warm-frontal  types  (extensive  area  of  weak  vertical 
vdocities  and  strong  mesoscale  ascent  over  a  more  limited  area,  respectively)  indicate 
that  in  the  first  case  seed  crystals  grow  by  vapor  deposition,  while  in  the  second  they 
grow  by  riming. 

Studies  of  precipitation  bands  from  other  regions  have  been  more  limited  in 
scope  than  the  (TYCILES  experiment.  Furthermcnre,  more  attention  was  given  to  the  role 
of  condititHial  symmetric  instability  (CSI)>  fust  proposed  by  Bennetts  and  Hoskins 
(1979)  and  investigated  further  by  Emanuel  (1983a,  b).  Reuter  and  Yau  (1990) 
examined  seven  cases  of  banded  precipitation  during  the  Canadian  Atlantic  Storms 
Program  and  concluded  that  slantwise  convection  may  have  been  a  factor  in  producing 
predpitation  bands  in  the  vidnity  of  Atlantic  low  pressure  areas  passing  south  of  Nova 
Scoda.  Wolfsberg  et  al.  (1986)  studied  precipitation  bands  on  the  order  of  50-100  km 
in  width  diat  occurred  in  a  New  England  winter  storm  and  found  that  they  paralleled  the 
thermal  wind,  characteristic  of  CSI.  However,  they  concluded  that  other  factors  also 
contributed  to  CSI  in  forming  and  maintaining  the  bands  including  frontogenetical 
forcing.  Mesoscale  bands  were  observed  in  the  coastal  front  during  a  GALE  case  study 
ccHiducted  by  Riordan  (1990).  These  bands  were  attributed  to  frcmtogentic  fordng, 
differential  heating,  and  differences  in  land-sea  friction.  Shields  et  al.  (1991)  identified 
three  types  of  snow  bands  in  a  winter  storm  over  east-central  Illinois  during  UNIWIPP. 
They  concluded  that  two  narrow  intense  bands  early  in  the  storm  were  the  result  of 
surface  ctmvergence  along  a  narrow  confluence  zone  associated  with  the  remnants  of  a 
weak  surface  warm  front  The  second  type  of  band  was  a  result  of  surface  convergence 
in  die  vicinity  of  an  involed  pressure  trough.  Multiple  parallel  bands  framing  in  the 
later  stages  of  the  observation  period  were  the  result  of  CSI  from  destabilization  caused 
by  frontogenetic  forcing.  Byrd  (1989)  performed  a  composite  analysis  for  27  winter 
overrunning  precipitation  events  over  the  southern  Great  Plains  and  found  CSI  could  be 
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a  fencing  mechanism  in  the  majority  of  the  cases  he  classified  as  banded  and  strongly 
banded.  Lindzen  and  Tung  (1976)  proposed  ducted  gravity  waves  have  a  roll  in  the 
devdt^ment  and  maintenance  of  precipitation  bands  observed  in  association  with 
convective  storms. 

Within  widespread  precipitation  regions  of  northeast  Colorado  winter  storms, 
heavier  snowfall  has  been  observed  to  coincide  with  mesoscale  bands  having  an 
assortment  of  orientations  (Rasmussen  et  al.  1992).  Lilly  (1981),  presented  radar 
observations  of  bands  parallel  to  the  mean  wind  vector.  Dunn  (1988)  examined  a  large 
band  associated  with  a  September  snowstorm  over  northeast  Colorado.  He  found  that 
(nogrtq)hic  influences  were  unimpentant  in  this  case  and  that  CSI  coupled  with 
ageostii^hic  jet  streak  and  frontogenetical  circulations  led  to  the  conditions  necessary 
for  fmnation  of  the  band.  Wesley  and  Pielke  (1990)  observed  banded  structure  in  two 
cases.  Similarities  in  the  data  between  their  cases  and  that  of  Wolfsberg  et  al.  (1986) 
(Le.  bands  parallel  to  the  thermal  wind  vector,  similar  sounding  structure,  and  radar 
observations)  lead  than  to  conclude  that  CSI  could  have  been  a  factor  in  band 
foimatitm.  However,  they  also  raised  the  question  of  the  roles  played  by  other 
mechanisms  such  as  low  level  convergence  zones,  mountain-induced  gravity  waves, 
and  lee  waves.  Rasmussen  et  al.  (1990)  documented  a  case  of  snow  bands  that  were 
not  parallel  to  the  thermal  wind  thus  failing  one  of  the  qualitative  tests  for  CSI.  Their 
modeling  studies  indicate  that  a  density  current  propagating  southward  with  a  cold  firont 
was  the  mechanism  responsible  for  band  growth. 

In  summary,  mesoscale  precipitation  bands  are  frequently  observed  in  winter 
storms.  A  large  body  of  evidence  suggests  that  CSI  may  play  significant  role  in  the 
femnation  and  nutintenance  in  these  bands.  However,  other  mechanisms  might  also  be 
iirqxntant  ctmtributors  to  snow  bands  over  northeast  Colorado.  They  include:  low  level 
convergence  zones,  low  level  barrier  jets  (Wesley  1991),  conditional  instability,  ducted 
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gravity  waves,  mountain  or  lee  waves,  upslope  induced  vertical  motions,  density 
currents,  and  upper  level  jet  streaks.  These  phenomena  will  be  evaluated  in  greater 
detail  for  their  role  in  band  formation  with  respect  to  this  case  in  Chapter  4. 

1,5  Radar  Data 

Radar  has  played  a  vital  role  in  the  observation  and  analysis  of  winter 
precipitation  bands,  including  the  studies  cited  above.  Several  other  studies  merit 
attention  here.  Szoke  (1991)  described  the  use  of  the  MHR  in  an  operational  setting 
where  snow  bands  were  observed  over  the  Denver  metropolitan  area  and  correspraided 
to  regions  of  heavier  precipitation.  Heckman  and  Dulong  (1989)  showed  there  was  a 
good  relationship  between  northeast  C>)lorado  snow  band  orientation  and  movement 
with  die  500  mb  height  and  wind  fields  respectively.  Sanders  (1986)  investigated  a 
major  New  England  snow  band  associated  with  rapid  development  of  an  offshore  low 
pressure  area.  He  found  frontogenetic  forcing  was  the  most  important  factor  in  band 
development 

Radar  data  analysis  techniques  vary  from  simple,  qualitative  examinatitm  of  Plan 
Position  Indicator  (PPI)  and  Range  Height  Indicator  (RHI)  displays  to  dual-Doppler 
studies  and  are  summarized  by  Battan  (1973)  and  Doviak  and  Zmic  (1984).  Most 
existing  studies  of  snow  bands  used  single-Doppler  analysis  techniques.  However, 
Heymsfield  (1979)  performed  a  dual-Doppler  analysis  of  precipitation  bands  parallel  to 
a  warm  finont  moving  over  Chicago  to  investigate  the  three  dimensional  structure  of  the 
bands.  The  bands  were  oriented  perpendicular  to  the  mid  level  flow,  contrary  to  bands 
examined  in  this  case  where  they  paralleled  the  upper  level  winds.  Dual-Doppler  study 
of  snow  bands  near  Champaign,  Ill  was  conducted  by  Bohne  (1979).  The  bands 
formed  in  the  rear  of  a  cyclone  and  in  general  were  perpendicular  to  the  flow,  but  nxrved 
in  the  direction  of  the  flow.  Kessinger  and  Lee  (1991)  describe  an  operational  dual- 
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Dq)pler  analysis  technique  used  during  WISP.  In  general,  the  literature  contains  little 
on  dual-Doppler  studies  of  snow  bands,  especially  over  northeast  Colorado.  Dual- 
Dq)pler  analysis  was  performed  on  one  volume  scan  in  this  case  to  examine  the 
kinematic  structure  of  small  scale  snow  bands. 

In  this  study  single-Doppler  analysis  provided  information  on  the  larger  scale 
organization  and  the  movement  of  the  storm  system,  major  snow  bands,  and 
precipitation  cores.  The  velocity-azimuth  (VAD)  display  developed  by  Browning  and 
Wexler  (1968)  has  been  used  to  investigate  the  kinematic  stmcture  of  horizontal  wind 
fields  in  winter  storms  (Lilly  1981,  Wolfsberg  etal.  1986,  Sanders  and  Bosart  1985b, 
Sassoi  et  al.  1990).  Here,  an  attempt  was  made  to  use  the  extended  velocity-azimuth 
di^lay  (EVAD)  technique  (Srivastava  et  al.  1986,  Matejka  and  Srivastava,  1991)  to 
retrieve  estimates  of  mesoscale  vertical  motion.  The  results  of  radar  observations  and 
analyses  are  presented  in  Chapter  3. 
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Figure  1.1:  WISP-91  Instrumentation  locator  map.  Dots,  crosses,  and  X’s  represent 
Si^ace  Airway  Observation  network  stations,  PAM  stations,  and  PROFS  stations, 
reflectively.  ^U-CHILL  and  Mile  High  radars  (indicated  by  CHL  and  MHR)  dual- 
Doippler  30  degree  crossing  angle  lobes  are  the  large  circles.  Dual-Doppler  analysis  is 
the  rectangular  region.  CLASS  soundings  were  t^en  at  PTL,  P03,  P12,  P14,  P16, 
P32,  P35,  and  near  PIO.  Wind  Profilers  and  microwave  radiometers  were  located  at 
PTL  and  DEN  and  PTL,  DEN  and  P32  respectively.  Terrain  contours  are  in  meters. 
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Figure  IJ2:  Terrain  height  contours  every  200  m  for  northeastern  Colorado.  The  2000, 
3000,  and  4000  meter  contours  are  indicated  by  bold  lines.  Major  Surface  Airway 
Obseivation  netwoiic  stations  are  also  shown. 


Chapter  2 

METEOROLOGICAL  ANALYSIS 

This  chapter  describes  the  synoptic  and  mesoscale  situation  pricn*  to  and  during 
the  snowstorm  of  19-20  January  1991.  In  general,  an  arctic  cold  front  moved  through 
eastern  Colcnado  leaving  anticyclonic  upslope  flow  conditions  favorable  for  the  initiation 
of  snow  (Boatman  and  Reinking  1984,  Reinking  and  Boatman  1987,  Toth  1991,  Wesley 
1991).  Observations  showed  snow  started  over  most  of  northeast  Colorado  around 
2100  UTC  19  January  and  ended  by  1200  UTC  20  January.  Total  snowfall  ranged 
between  less  than  2.5  cm  (1.0  in)  to  more  than  IS  cm  (6.0  in)  with  liquid  water 
equivalents  between  .03  cm  (0.01  in)  and  .75  cm  (0.30  in).  Figures  2.1  and  2.2  show 
aiowfall  and  liquid  water  equivalent  contours  respectively,  obtained  from  the  special 
WISP-91  snow  spotter  observer  network  sites  (Wesley  1991),  and  National  Weather 
Service,  Federal  Aviation  Administration,  and  Department  of  Defense  weather  stations. 
Note  the  three  distinct  snowfall  maxima;  the  tongue  extending  from  the  Fort  Collins  area 
across  Weld  County  and  into  M^gan  County,  the  region  north  of  Denver,  and  area  over 
the  eastern  Palmer  Divide.  A  radar  loop  of  Mile  High  radar  scans  and  surface 
observatitMis  suggest  that  topographic  influences  on  the  low  level  flow  regime  played  a 
significant  role  in  the  initiaticm  and  maintenance  of  heavier  snowfall  over  the  Denver 
region  and  the  maximum  near  Ftnt  Collins.  For  the  other  region,  the  upper  level  flow 
regime  played  a  nxxe  important  role,  particularly  over  the  Palmer  Divide  where  a  well 
defined  snow  band  was  present.  More  discussion  of  these  two  flow  regimes  and  radar 
observations  will  be  given  in  chapter  3. 
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Unlike  the  arctic  outbreak  case  examined  by  Wesley  (1991)  which  occurred  over 
5  days,  no  nxne  than  18  hours  elapsed  from  frontal  passage  until  the  end  of  the  upslope 
generated  snow.  Also,  there  were  other  departures  from  the  anticyclonic  upslope 
situation  described  by  Boatman  and  Reinking  (1984).  These  differences  are  addressed 
below  and  in  Chapter  3. 

2.1.  Upper  Air 

A  strong,  vigorous  shrHtwave  in  the  upper  atnnosphere  and  associated  arctic 
surface  cold  front  moved  southeast  out  of  Canada  on  19  January  passing  to  the  east  of 
Colorado  by  1200  UTC  20  January.  The  upper  level  shortwave  trough  had  its  origins 
over  the  North  Pacific  but  came  in  phase  with  a  polar  low  pressure  area  west  of  northern 
Greenland  that  had  a  southeast  extension  over  Hudson  Bay.  National  Meteorological 
Center  (NMC)  500  mb  analyses  (Figs.  2.3a-d)  for  19  and  20  January  revealed  the 
evtdution  and  movonent  of  the  upper  level  trough.  During  this  period,  the  Denver  500 
mb  height  decreased  180  meters,  while  the  l(X)0-500  mb  thickness  decreased  from  over 
5520  to  less  than  5280  meters.  Thermal  wind  analyses  were  constructed  fen*  various 
ccHnbinations  of  sounding  levels.  Figure  2.4  is  representative  of  these  analyses  and 
indicated  a  southwest  to  northeast  orientation  to  the  800-4(X)  mb  thermal  wind.  Radar 
obsmved  precipitation  bands  (described  in  Chapter  3),  particularly  above  the  boundary 
layer,  paralleled  the  thermal  wind,  characteristic  of  conditional  symmetric  instability 
(CSI)  (Sanders  and  Bosart  1985). 

NMC  vorticity  analyses  fOT  the  period  (not  shown)  indicated  a  well  defined 
vorticity  maximum  associated  with  the  shortwave  trough  moved  across  the  state  from  the 
ntHthwest  comer  to  the  southeast  comer  between  (X)(X)  and  12(X)  UTC  20  January. 
Northeast  Colorado  was  in  an  area  of  positive  vorticity  advection  by  0000  UTC  20 
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January.  The  500  mb  temperatures  over  Denver  decreased  12'C  from  -18  to  -30*C, 
during  the  period  which  indicated  strong  cooling  aloft. 

Winds  increased  and  backed  with  dme  as  the  jet  stream  migrated  over  the  regicm 
between  0000  UTC  19  January  and  1200  UTC  20  January  and  was  reflected  in  the  NMC 
300  mb  arudysis  for  the  period  (Figs  2.5a-d).  Two  branches  of  the  jet  stream  were 
visible  at  (XXX)  UTC  19  January  1991  with  split  flow  around  the  low  pressure  area  over 
Texas.  A  strong  jet  streak  at  1200  UTC  19  January  1991  associated  with  the  shortwave 
over  Idaho  and  western  Washington  weakened  and  shifted  eastward  by  0000  UTC  20 
January  1991.  Meanwhile,  the  southern  Iranch  of  the  jet  stream  shifted  northward  and 
merged  with  the  northern  portion  of  the  split  flow  as  the  500  mb  low  pressure  area 
weakened  over  Oklahoma.  At  0000  UTC,  during  the  height  of  the  storm,  northeast 
Colorado  was  not  in  either  the  classic  left-exit  or  right-entrance  quadrant  region  of  a  well 
defined  jet  streak,  the  favored  region  for  heavy  snow  (Dunn  1988).  Examination  of  the 
higher  temporal  and  spatial  resolution  offered  by  the  3(X)  mb  CLASS  sounding  plots 
(Hg.  2.6a-d)  supports  this  conclusion.  Reinking  and  Boatman  (1986)  suggest  that 
precipitaticxi  genoated  by  terrain  induced  upward  vertical  motions  of  the  arctic  airmass 
can  occur,  despite  subsidence  associated  with  the  left-entrance  region  of  jet  streaks. 
Furthermore,  tenqrerature  advection  can  modify  the  favored  locations  of  jet  streak 
induced  midtropospheric  vertical  modon  (Keyser  and  Shapiro  1986). 

Hme  series  of  the  Platteville  and  Denver  windprofilers  (Figs.  2.7a-b)  depict  the 
vertical  structure  of  horizontal  winds  at  every  hour  between  1200  UTC  19  January  and 
0600  UTC  20  January  1991.  Several  features  can  be  seen: 

1.  A  wind  maximum  was  present  in  the  lower  levels  of  the  Platteville  profiler  in 
the  early  observations  prior  to  passage  of  the  surface  cold  front,  indicating  the  presence 
of  moderate  downslope  winds.  This  was  confirmed  by  surface  observations  at  Fort 
Collins  where  a  gust  of  25  ms*^  from  the  west  was  recorded  at  1255  UTC. 
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2.  An  upper  level  wind  maximum  crossed  the  region  between  2000  UTC  19 
January  and  0300  UTC  20  January. 

3.  The  20  minute  resolution  Denver  profiler  plot  (Figure  2.7b)  showed 
deepening  and  backing  of  low  level  northeasterlies  associated  with  the  surface  arctic 
airmass.  This  was  also  visible  in  the  CLASS  time  series  (Figures  2.14a-d). 

4.  ^^nds  increased  with  height  and  had  very  little  directional  shear  above  the 
boundary  layer.  Both  these  condidons  are  qualitative  indicators  of  conditional  symmetric 
instability  (Snook  1992;  Wolfsberg  etal.  1986). 

Examination  of  the  700  mb  moisture  field  (see  Figs  2.8a-c)  showed  a  significant 
portion  of  the  storm  moisture  was  carried  along  with  the  upper  level  trough.  I>ewpoint 
depressions  at  0000  UTC  19  January  (Fig.  2.8a)  over  most  of  the  western  U.  S.  were 
greater  than  12*C  while  those  in  the  vicinity  of  the  upper  level  trough,  while  variable 
over  space,  were  less  than  10*C.  These  lower  dewpoint  depressions  associated  with  the 
trough  could  be  traced  back  to  even  earlier  NMC  analyses  (not  shown).  Hgs.  2.8b-c 
show  the  moisture  at  700  mb  moved  with  the  trough.  Dewpoint  depressicHis  were 
slightly  Iowa*  at  500  mb  and  were  probably  the  result  of  moisture  advected  from  the 
upper  level  low  pressure  area  that  moved  over  C^ifomia  (see  Figs  2.3a-d).  Minimum 
mid-level  dewpoint  depressions  were  actually  above  the  500  mb  mandatory  reporting 
level  on  the  Denver  and  Grand  Junction  1200  UTC  19  January  soundings  (Figs  2.9a-b ). 
This  elevated  moisture  regitxi  was  also  visible  on  the  1800  UTC  19  January  Berthoud 
CLASS  sounding  (Hg  2.10a).  Boatman  and  Reinking  (1984)  point  out  the  need  for 
relatively  moist  Pacific  air  to  cross  the  Rocky  Mountains  and  Colorado  Front  Range  for 
cloud  enhancement  (and  greater  snowfall)  above  the  shallow  arctic  airmass  over  the 
Eastern  Rains.  In  this  case,  the  cross  mountain  rrxristure  requirement  was  actually 
satisfied  by  complex  interaction  of  the  two  distinct  sources  described  above. 
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By  0000  UTC  20  January  the  entire  airmass  had  approached  saturation  up  to  the 
tropopause,  which  lowered  nearly  70  mb  in  6  hours  at  Berthoud  (see  Figs.  2.10a-b). 
Similar  reductions  in  the  tropopause  height  occurred  at  all  ClASS  stations.  The  (XXX) 
UTC  Akron  and  Elbert  CLASS  soundings  (Figs,  2.10c-d)  resembled  the  Berthoud 
sounding  with  the  exception  of  a  shallower  cold  air  depth  at  Elbert  Comparison  of  the 
pressure  at  the  top  of  the  arctic  airmass  inversion  recorded  in  the  Berthoud  and  Akron 
soundings  suggested  observational  verificatitxi  of  Wesley's  (1991)  description  of  a 
vertical  bulge  in  the  western  edge  of  the  arctic  air  mass.  A  greater  easterly  wind 
component  in  the  boundary  layer  would  have  made  this  feature  more  pronounced  and 
possibly  enhanced  precipitation  near  the  Frcxit  Range  (see  Wesley  1991). 

Auer  and  White  (1982)  found  that  a  6(X)  mb  temperature  between  -12  and  -16*C 
was  a  common  characteristic  in  75  heavy  snow  events.  They  deduced  heavy  snow  will 
occur  wdien  the  temperature  regime  of  maximum  crystal  growth  by  vapor  deposition  in  a 
water  saturated  envircmment  (between  -12  and  -16*Q  coincides  with  maximum  vertical 
velocities  associated  with  the  level  of  nondivergence  (approximately  600  mb).  All  0000 
UTC  20  January  CLASS  soundings  had  600  mb  temperatures  between  -12  and  -16'C 
which  agreed  with  the  findings  of  Auer  and  White.  While  heavy  snow  did  occur  in  this 
case,  the  distribution  was  uneven  (see  Fig  2.1)  reflecting  the  role  of  other  processes  (e.g. 
snow  bands,  terrain  effects,  etc.)  in  local  snowfall  amounts. 

2.2  Vertical  Cross  Sections 

Cross  sections  in  space  and  time  of  various  data  and  parameters  were  constructed 
to  gain  understanding  of  the  thermodynamic  structure  of  the  atmosphere  during  this 
jnecipitaticm  event  Figure  2.1 1  shows  the  two  (X)00  UTC  20  January  vertical  cross 
sectimi  transects  discussed  in  this  section  as  well  as  the  low  level  winds  at  the  stations  in 
Northeast  Colorado.  The  arctic  airmass  was  visible  as  the  tighter  packing  of  isentropes 
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in  the  lower  levels  of  Fig.  2.12a ,  and  between  DEN  and  LBF  (see  Appendix  A  for  a 
listing  of  statitm  identifiers,  names,  and  location  information)  in  Fig.  2. 12b.  The  cross 
section  planes  in  Figures  2. 12a  and  2.12b  were  oriented  approximately  perpendicular 
and  within  30  degrees  of  parallel  to  the  low  level  winds  respectively.  In  the  lower  levels 
these  cross  sections  were  similar  to  those  of  Boatman  and  Reinldng  (1984),  with  the 
exception  that  isentropes  sloped  downwards  between  BTD  and  DEN.  Additionally, 
above  die  arctic  airmass,  a  tendency  towards  lower  potential  stability  values  reflected  a 
more  vigorous  shortwave  aloft  than  observed  by  Boatman  and  Reinldng  (1984).  Figure 
2.12a  also  shows  that  the  arctic  airmass  had  just  passed  Dodge  City  as  reflected  in  the 
position  of  the  surface  front  in  Fig.  2.1Sb. 

Examination  of  the  time  series  plots  of  equivalent  potential  temperature  at  the 
CLASS  stations  Akron  (AKR),  Wiggins  (WIG),  Berthoud  (BTD),  and  Flaglo-  (FLA) 
(Rgs.  2.13a-d)  showed  conditional  instability  near  the  surface  at  BTD  and  FLA,  and  a 
trend  towards  instalnlity  above  600  mb  in  all  cross  sections.  Increasing  relative 
humidities  in  the  cross  sections  were  more  the  result  of  cold  air  advection  rather  than 
significant  moisture  advection.  The  CLASS  sounding  data  showed  only  slight  increases 
in  the  dewpoint  ten:q)etature  profiles,  whereas  there  was  significant  cooling  through  the 
depth  of  the  atmosphere  with  the  passage  of  the  surface  cold  front  and  subsequent  arrival 
of  arctic  air  in  die  low  levels,  and  cold  air  advection  associated  with  the  upper  level 
trough  above  the  arctic  airmass. 

Whiteman  (1973)  states  that  over  the  high  plains,  relative  humidities  near  85% 
define  regions  of  cloudiness  and  anticyclonic  steams  are  often  characterized  by  abrupt 
decreases  in  humidity  (increased  dewpoint  depressions)  at  cloud  te^  as  opposed  to  the 
more  gradual  dewpoint  depression  increases  in  height  associated  with  deqi  cyclonic 
stonns.  Dewpoint  depressions  in  the  CLASS  soundings  (Figs.  2.10b-d)  increased 
(hamatically  at  the  tropopause  above  the  400  mb  level  indicating  a  distinct  cloud  top. 
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lending  support  to  the  conclusion  that  this  storm  can  be  classified  as  anticycltHiic 
(Rasmussen  et  al.  1992). 

It  has  been  proposed  that  given  the  proper  orientation  of  the  upper  level  flow  with 
reflect  to  die  frontal  boundary  aloft,  upward  vertical  motion  may  occur  in  the  flow  just 
above  the  boundary  up  to  several  hours  after  passage  of  a  shallow  sloped  surface  cold 
fxtmt  (personal  communication.  Dr.  R.  Rasmussen,  1992).  This  might  help  explain  the 
lag  in  cloudiness  and  precipitation  after  passage  of  the  surface  cold  front  often  observed 
with  winter  storms  over  northeast  Colorado.  A  rough  approximation  of  the  frontal  slope 
in  this  case  was  obtained  by  taking  an  average  of  the  vertical  displacements  of  equivalent 
potential  tonperature  surfaces  near  the  top  of  the  arctic  airmass  from  the  New  Raymer 
and  Elbert  (X)00  UTC  20  January  CLASS  sounding  data.  The  New  Raymer  Elbert  axis 
is  ai^noximately  parallel  to  the  directicm  of  movement  experienced  by  the  low  level  arctic 
airmass.  The  average  vertical  displacement  was  500  m  across  a  hcHizontal  distance  of 
165  km.  This  yields  a  slope  of  approximately  3.03  x  10*3  or  an  angle  of  0.17*  for  the 
arctic  airmass  ^ich  had  an  approximate  depth  of  1500  m  and  top  between  3^00-3900  m 
MSL  over  northeast  Colorado.  If  the  winds  above  the  top  of  the  boundary  were  from 
260  degrees  at  15  ms'l  (see  Fig.  2.10b-d),  and  the  axis  of  the  frontal  slope  was  oriented 
030-210*  (increasing  in  height  towards  030*),  then  15  ms*l  x  cos  50*,  or  9.6  ms*^  will 
represoit  the  component  of  the  wind  forced  upward  along  the  frontal  surface.  If  the 
frontal  surface  was  moving  at  a  speed  of  17  ms'^  (see  section  2.3),  then  the  vertical 
vdoci^  induced  the  upper  level  wind  flowing  up  the  arctic  airmass  boundary  will  be 
(9.6  +  17)  X  sin  0.17*,  or  0.08  ms*l  (8  cm  s*l).  This  is  probably  close  to  a  theoretical 
maximum  as  it  is  likely  that  the  airmass  had  slowed  to  less  than  17  ms'l.  As  described  in 
secticm  2.3,  the  surface  front  was  well  through  nonheast  Colorado  by  0000  UTC.  Thus, 
the  data  suggest  duit  conditions  were  favorable  for  weak  vertical  ascent  of  the  upper  level 
winds  which  may  help  explain  the  lag  in  the  onset  of  cloudiness  and  precipitation. 
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2.3  Surface 


In  response  to  the  advance  of  the  upper  level  trough  of  low  pressure  described 
above,  a  surface  low  pressure  area  and  associated  frontal  system  developed  in  Alberta. 
The  low  pressure  area  skirted  the  northern  tier  states  while  the  cold  front  moved  rapidly 
south  out  of  Canada,  through  the  northern  plains  states,  and  into  Northeast  Colorado 
over  an  18  hour  period.  This  represents  an  average  speed  of  nearly  60  km  hr'^  (or  ~17 
ms*^).  According  to  Reinking  and  Boatman  (1987)  this  sequence  of  events  is  typical  fn- 
anticylconic  upslope  conditions  over  N(»iheast  Colorado.  However,  at  the  same  time  a 
trough  of  low  pressure  over  eastern  Colorado  and  the  subsequent  development  of  a 
surface  low  pressure  area  over  southeast  Colorado  in  response  to  the  upper  level  trough, 
strengthened  the  post  frcmtal  pressure  gradient  over  the  eastern  plains.  This  aided  in  the 
development  of  upslope  flow  rather  than  having  to  rely  on  post  frontal  anticyclraiic  flow 
alone.  Figs.  2.15a>b  show  the  surface  synoptic  situation  as  analyzed  by  the  National 
Meteorological  Center  (NMQ  at  1200  UTC  19  January  and  0000  UTC  20  January 
1991.  They  illustrate  the  speed  with  which  the  front  moved  and  development  of  the 
southeast  Qdorado  low  pressure  area.  Prior  to  frontal  passage  northern  portions  of  the 
Colorado  fnmt  range  experienced  downslope  wind  conditions  with  gusts  to  25  ms'^  in 
respcmse  to  die  surface  pressure  gradient  between  high  pressure  over  Utah  and  a  lee  side 
trough  over  eastern  Colorado  and  western  Kansas. 

Hme  cross  sections  were  plotted  for  all  PAM  and  PROFS  stations  used  in  die 
experiment.  The  Greeley  cross  section  (Fig.  2.15)  is  typical  of  conditions  experienced  at 
most  stations  befne,  during,  and  after  frontal  passage.  In  this  case  there  was  no 
dramatic  rise  in  pressure  and  only  a  moderate  temperature  decrease,  which  omtrasts  with 
the  cases  desoibed  by  Marwitz  and  Day  (1991),  Wesley  (1991),  and  the  dramatic 
exan^le  of  Shapiro  (1984)  where  rapid  temperature  decreases  occurred  shordy  after 
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finontal  passage.  A  similar  gradual  surface  temperature  decrease  occurred  in  the  case 
described  by  Dunn  (1987).  However,  in  that  case  surface  winds  had  a  much  greater 
easterly  conqmnent  which  lead  to  cold  air  damming  and  the  development  of  a  strong 
convergence  zone,  or  meso-front,  some  distance  away  from  the  rapidly  rising  terrain  of 
the  Front  Range.  The  greatest  temperature  falls  in  the  case  studied  here  took  place  over 
an  hour  after  the  1900  UTC  frontal  passage,  as  indicated  by  a  wind  shift  and  increased 
windspeeds.  Snow  began  around  the  time  that  the  temperature  -  dewpoint  spread 
decreased  significantly.  The  Greeley  Airport  weather  observer  took  a  special  observation 
for  snow  at  2124  UTC  which  is  shortly  after  the  large  dewpoint  rise.  At  2146  UTC  the 
Denver  Airport  National  Weather  Service  Forecast  Office  reported  snow.  It  continued 
there  until  1041  UTC  20  January  with  the  heaviest  amounts  falling  between  2200  UTC 
19  January  and  0600  UTC  20  January.  Radar  data  showed  that  the  onset  of  snow  was 
n^nd,  occurring  ovct  most  of  Northeast  Colorado  within  a  shOTt  time  around  2100  UTC, 
exo^t  fix  a  region  in  the  iKMtheast  pxtion  of  the  radar  coverage  area  v  'bere  snow  began 
shortly  after  frontal  passage.  Surface  data  it^cated  that  dewpoints  decreased  at  a  slower 
rate  dian  the  tanperature  after  frontal  passage  which  suggests  that  cold  air  advectitxi 
rather  than  influx  of  moisture  laden  air  from  some  other  region  was  the  cause  for 
increasing  low  level  humidities  after  frontal  passage.  This  cooling  process,  in 
ccmjunction  with  similar  cooling  described  in  section  2.1,  todc  several  hours  which 
offers  a  partial  explanation  for  the  lag  in  precipitation  after  passage  of  the  surface  front.. 

Unlike  the  cases  described  by  Wesley  and  Pielke  (1990),  Dunn  (1987),  and 
Wesley  (1991),  there  was  no  evidence  to  support  an  argument  for  cold  air  damming.  No 
distinct  cold  pool  along  the  foothills  and  adjacent  plains  was  present  in  the  isotherm 
analyses  (sec  Rg  2.16  for  example).  However,  colder  temperatures  occurred  over  the 
higher,  mountainous  terrain  to  the  west  and  well  into  the  arctic  airmass  to  the  nenth. 

Also,  time  was  no  indication  of  a  barrier  jet,  common  in  extreme  cases  of  cold  air 
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damming  (Wesley  and  Pielke  1990,  Wesley  1991),  in  the  wind  profiler  (see  Figs.  2.7a- 
b)  or  surface  wind  data. 

Figure  2.17  shows  a  typical  streamline  and  isotach  pattern  over  the  region  during 
the  storm.  Note  the  only  evidence  of  a  confluence  zone,  and  thus  meso-fiont  as 
described  by  Wesley  and  Pielke  (1990),  Wesley  (1991),  Boatman  and  Reinking  (1984), 
and  Dunn  (1987)  was  over  the  area  south  of  Denver  and  near  the  western  portions  of  the 
Palma*  Divide.  Pure  upslope  flow  was  a  significant  factor  in  low-level  upward  vertical 
motion,  especially  over  the  Palmer  Divide  (and  a  small  region  near  Fort  Collins). 
Assuming  there  is  an  average  700  m  rise  in  terrain  over  a  distance  of  approximately  125 
km  fixm  the  Platte  River  Valley  east  of  Greeley  to  the  ridge  of  the  Palmer  Divide  near 
PAM  station  P31,  upslope  winds  with  an  upslope  component  of  7  ms"^  (or  ~25  km 
hr*)  would  lead  to  a  vertical  velocity  of  nearly  4  cms**.  Therefore,  upslope  flow  and 
low-level  convergence  may  partially  explain  the  higher  precipitation  over  the  region  south 
of  Denver.  Also,  easterly  flow  in  the  vicinity  of  Fort  Collins  provided  strong  upslope 
flow  and  probably  played  an  important  role  in  causing  heavier  snowfall  in  that  region. 
Other  faaors  (e.g.  seeder-feeder  mechanism,  CSI,  etc.)  that  may  have  also  been 
responsible  for  the  fexmation  and  maintenance  of  snow  bands  observed  by  radar  and 
respcmsible  for  the  uneven  distribution  of  precipitation  will  be  exnlored  in  later  chapters. 

2.4  Summary 

Anticyclonic  upslope  flow  in  the  wake  of  an  upper  level  trough  of  low  pressure 
and  associated  frcxital  system,  and  enhanced  by  the  development  of  a  surface  low 
pressure  area  in  southeastern  Colorado,  brought  5-10  cm  of  snow  (over  15  cm  in  few 
locations)  to  Northeast  Colorado.  The  system  appeared  to  have  two  distinct  moisture 
sources  in  the  upper  atmosphere,  one  near  700  mb  that  traveled  with  the  trough  and  the 
otiier  having  origins  in  an  upper  level  low  pressure  area  over  C^ifomia.  Cold  air 


24 


H^iinming  and  a  bairier  jet  were  not  present  during  this  storm.  There  was  some 
confluence  in  the  wind  south  of  Denver,  as  well  as  upslope  conditions  near  Fort  Collins, 
but  this  does  not  fully  explain  the  precipitation  distribution  pattern  or  the  highly  banded 
structure  of  the  precipitation  observed  by  radar. 
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Figure  2.1;  Total  snowfall  distribution  (in)  over  Northeast  Colorado  19-20  January 
1^1.  Terrain  height  contours  are  every  5(X)  m. 


Figure  2.2:  Total  liquid  water  equivalent  (in)  over  Northeast  Colorado  19-20  January 
1^1.  Terrain  height  contours  are  every  500  m. 
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D155  ..  SOOMB  ANALYSIS  HE  I 


Figure  2.4: 0000  UTC  800-400  mb  thermal  wind  streamlines  (solid  lines)  and  isotachs 
(dadied  lines)  in  ms*^  Line  represents  axis  along  which  CSI  profile  was  computed  in 
Fig.  4.3. 
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analysis  heights/ I  sot  ACHS  00/  SUN  20  JAN  9\ 


Rgure  2.6:  QASS  300  mb  station  plots  over  Ncothcast  Colorado  for  a)  1800  UTC  19 
January  1991,  b)  2100  UTC,  c)  0000  UTC  20  January,  and  d)  0^  UTC.  From  upper 
left  comer  of  plot  working  clockwise,  plots  show  temperature  in  *C,  height  in  dm, 
station  identimr,  luid  dewpdnt  in  *C.  Pennants,  long  barbs,  and  short  barbs  on  wind 
riiafts  rqnesent  25  ms*',  5  ms*',  and  2.5  ms*'  respectively. 
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lALYSIS  HEIGHTS/TEMPERATURE  OOZ  SAT  19  JAN  91 
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D27B  ..  700MB  ANALYSIS  HEIGHTS/TEMPERATURE  12Z  SAT  19  JAN  91 


54 


IS  Z00 


a) 


910119/1200  72^:69  OCN 


910119/1200  72<I76  GJT 


Rgure  2.9: 1200  JJTC  19  January  1991  rawinsonde  sounding  for  a)  Denver  and  b) 

sohd  heavy  dashed  line,  horizontal  solid  Une.  straight 
sltgMy  curv^  hnes  and  heavily  curved  lines  represent  temperature  trace 
dCTijNMnt  tracysobars  (mb),  isothenns  (  C),  dry  adiabats,  and  moist  adiabS 
respectively.  TTic  same  wind  baib  notation  as  Fig.  2.7a  is  used  here 


56 


and  d)  Elbert 


Hg^  2.1 1:  Cnoss  section  transects  and  representative  low  level  winds.  Lcme  barbs  on 
wind  shafts  represent  5  ms-*,  and  short  barbs  2.5  ms'l. 


59 


60 


Figure  2.12:  (a)  0000  UTC  20  Januaiy  1991  Lander,  WY  to  Dodge  City,  KS  cross  section.  Solid 
lines  are  isentropes  (*K)  and  dashed  lines  are  relative  humidities  (percent).  Pennants,  long  barbs, 
and  short  barbs  on  wind  shafts  represent  50  ms’*.  10  ms*^  and  5  ms'*  respectively. 
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Figure  2.12:  (b)  Same  as  (a)  except  between  Grand  Junction,  CO.  and  North  Platte, 
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OOZ  S(K  20.  JAN  1991 


Figure  2.16: 0000  UTC  20  January  surface  isotherm  analysis  over  Northeast  Colorado. 
Heavy  lines  are  isothenns  fC)  and  thin  lines  are  500  m  terrain  height  concurs. 


Figure  2.17: 0100  UTC  20  January  surface  streamline  and  isotach  analysis  over 
NordMcast  Colorado.  Lines  with  arrows  are  streamlines,  heavy  dashed  Unes  are  speeds 
in  ms'l,  and  diin  lines  are  500  m  terrain  height  contours. 
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Chapter  3 

RADAR  OBSERVATIONS  AND  ANALYSIS 

In  this  chapter  attention  is  given  to  description  of  observations  and  analysis  of 
single-  and  dual-Dc^pler  radar  data.  Previous  studies  on  precipitation  bands  in  winter 
storms  cited  in  Chapter  1  made  extensive  use  of  radar  data  to  locate,  measure  intensities, 
and  evaluate  kinematic  properties  of  the  bands.  Similarly,  use  of  WISP-91  radar  data 
aided  in  the  observation  and  analysis  of  snow  bands  occurring  during  the  19-20  January 
storm.  Organized  patterns  in  PPl  displays  of  equivalent  reflectivity  (Ze)  were  used  to 
identify  the  bands,  their  orientation,  and  their  structure.  A  brief  description  of  data 
collectitm  and  processing  techniques  is  given  in  section  3.1.  Single-Doppler  radar 
observations  documenting  bands  of  several  scales,  orientations  both  parallel  and 
perpendicular  to  the  upper  level  flow,  and  precipitation  cores  are  presented  in  section 
3.2.  Estimates  of  mesoscale  vertical  motion  (discussed  in  section  3.3)  were  made  using 
die  extended  vertical-azimuth  display  (EVAD)  technique  (Srivastava  et  al.  1986,  Matejka 
and  Srivastava  1991).  Finally,  section  3.4  describes  a  dual-Doppler  analysis  performed 
at  2327  UTC  19  January  to  describe  a  particular  snow  band  in  greater  detail.  Appendix 
C  lists  the  radar  characteristics  for  the  Colorado  State  University  (CSU)  CHILL  and 
National  Oceanic  and  Atmospheric  Administration  (NOAA)/NCAR  Mile  High  (MHR) 
radars. 

3.1  Data  Collection  and  Processing 


71 


Data  collection  at  the  CSU-CHILL  radar  began  at  181 1  UTC  19  January  and 
lasted  until  0030  UTC  (transmitter  arcing  forced  an  end  to  operations  at  this  dme).  The 
MHR  collected  data  between  1719  UTC  19  January  and  1036  UTC  20  January.  During 
operations,  scientists  at  CSU-dHILL  monitored  the  storm  progress  and  communicated 
with  the  WISP  operations  center  in  Boulder,  CO  by  radio.  Decisions  on  scanning 
strategy  were  made  on  the  basis  of  on-site  scientist  inputs,  support  for  aircraft 
operations,  and  storm  type.  For  this  case  the  CSU-CHILL  alternated  between  360* 
volume  scans  consisting  of  24  elevation  angles,  with  the  highest  being  22.5*,  and  RHI 
scans.  In  an  effort  to  keep  the  size  of  radar  data  sets  more  manageable,  data  were  not 
collected  for  the  portion  of  any  beam  that  was  above  6  km  AGL  (which  was  above 
observed  echo  tops  for  this  event).  The  scanning  strategy  of  the  MHR  consisted  of  21 
elevation  angle  scans,  with  the  highest  being  17  degrees.  All  CSU-CHILL  volume 
scans  began  at  the  same  time  as  MHR  volume  scans. 

Field  data  from  CSU-CHILL  radar  were  recorded  on  nine  track  tapes  for 
subsequent  transfer  to  a  universal  format  (NCAR  1991)  Exabyte®  tape  for  processing 
on  Sun®  workstations.  MHR  data  were  archived  on  the  NC^R  Mass  Storage  System 
(MSS)  after  con^letion  of  the  experiment.  Radar  volurtKS  of  interest  were  retrieved 
from  MSS  for  local  processing  by  file  transfer  protocol  (ftp)  over  the  high  speed 
INTERNET  data  netwOTk.  Volume  scans  between  19(X)  and  2356  UTC  19  January  and 
2104  UTC  19  January  to  0053  UTC  20  January  were  processed  for  CSU-CHILL  and 
MHR  radars  respectively.  Appendix  B  summarizes  the  major  processing  steps 
necessary  to  edit  folded  velocities  and  noisy  data,  conven  spherical  coordinate  data  to 
Cartesian  coordinates,  and  display  the  results  as  constant  altitude  plan  posidcm  indicator 
(CAPPI)  di^lays  or  vertical  cross  sections,  or  use  the  (Dartesian  data  to  generate  dual- 
Dofqiler  syndieses.  Data  were  used  in  radar  polar  coordinates  for  EVAD  analysis. 
Regions  of  Dopplor  velocity  folding  were  small  due  to  low  target  velocities  with  respect 


72 


to  the  Nyquist  velocities  (28.67  ms'l  and  24.75  ms'l  for  the  CSU-CHILL  radar  and 
MHR  respectively).  Data  contaminated  by  ground  clutter  near  the  radar  or  returns  from 
the  Rocky  Mountains,  approximately  40-60  km  southwest  to  northwest  of  the  radars 
(see  Rg.  1.1)  were  deleted  using  the  Research  Data  Support  System  (RDSS)  (Mohr  et 
al.  1986).  Data  points  with  radial  velocities  ±  2  ms'l  and  reflectivities  exceeding  35 
dBZ  (precipitation  echoes  were  generally  less  than  35  dBZ)  within  80  km  of  the  radar, 
were  deleted,  though  care  was  taken  to  avoid  eliminating  good  data.  After  the  first  few 
volumes,  regions  that  corresponded  to  ground  clutter  or  mountain  returns  were  known 
with  sufficient  accuracy  that  only  the  bad  data  were  removed. 

3.2  Single>Doppler  Observations 

Passage  of  the  cold  front  between  1800  and  2100  UTC  19  .Tanuary  over 
nortiieast  Colorado  signaled  the  onset  of  north-northeast  flow  and  upslope  conditions 
over  the  region,  particularly  over  the  northern  portion  of  the  Palmer  Divide.  Echoes 
developed  rapidly  between  1900  and  2100  UTC.  After  2100  UTC  widespread 
jnecipitation  goterated  neariy  solid  radar  echoes,  particularly  below  3  km  MSL 
(^rproximately  1600  m  and  14(X)  m  AGL  at  CSU-CHILL  and  MHR  respectively)  and 
over  the  northern  portions  of  the  radar  coverage  area.  Above  3  km,  echoes  took  on  a 
distinct  banded  structure  that  paralleled  the  8(X)-400  mb  thermal  wind  field.  This 
organization  to  the  echo  pattern,  with  the  exception  of  variations  due  to  movement  and 
intensity  changes,  persisted  through  at  least  0100  UTC  when  radar  data  processing  was 
discontinued.  Figures  3.  la-d  show  a  series  of  MHR  2249  UTC  CAPPI  reflectivity 
plots,  which  charac^ze  the  echo  cvganization  during  the  observation  period,  ftn*  the 
levels  of  2.5, 3.5, 4.5,  and  5.5  km  (MSL).  Snow  was  widespread  over  the  northern 
portions  of  the  region  at  the  lower  levels  with  no  obvious  echo  structure.  At  higher 
levels  and  particularly  to  the  south  of  the  radar,  echoes  exhibit  a  highly-banded 
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structure.  Bands  were  (»iented  south-southwest  to  north-northeast,  and  paralleled  the 
upper  level  winds.  This  is  similar  to  the  cases  described  by  Wolfsberg  et  cd.  (1986), 
Shields  et  cd.  (1991),  Seltzer  et  al.  (1985)  and  others. 

The  radial  velocity  field  corresponding  to  Figs.  3.1a-d  is  shown  in  Figs  3.2a-d. 
The  2.5  km  (MSL)  CAPPI  shows  well  defined  north-northeast  flow.  Soundings  and 
wind  profiler  data  indicated  the  depth  of  this  layer  to  be  on  the  order  of  1.5  to  2  km. 
There  was  little  change  in  this  depth  over  the  radar  observational  period.  Winds  backed 
r^idly  over  a  shallow  shear  layer  near  4  km  (MSL)  and  became  west-southwesterly  and 
increased  in  speed  with  height  above  4  km  (MSL).  There  was  little  change  in  wind 
direction  above  the  shear  layer. 

After  2330  UTC  19  January  low  level  reflectivities  developed  and  increased  in 
area  over  a  region  north  of  the  western  portions  of  the  Palmer  Divide.  The  growth  and 
intensificaticm  of  low  level  echoes  is  illustrated  by  the  reflectivity  field  changes  at  2.5  km 
observed  between  2346  UTC  19  January  and  0044  UTC  20  January  (Figs.  3.3a-b).  A 
video  loop  of  the  MHR  1.6  degree  elevation  scans  for  the  period  between  19(X)  UTC  19 
January  and  0300  UTC  20  January  showed  the  region  grew  slowly  in  areal  extent,  and 
moved  south-southwest  coincident  with  the  low  level  north-northeast  winds.  Periods  of 
moderate  snow  at  Jefferson  County  Airport  (BJC),  Ontennial  Airport  (APA),  and 
Denver  (DEN)  corresponded  to  passage  of  this  low  level  enhanced  reflectivity  region. 
Centomial  Airport  (APA),  the  station  closest  to  the  northern  slope  of  the  Palmer  Divide, 
reported  naxlerate  snow  between  0145  and  06(X)  UTC  20  January  while  Jefferson 
County  Airport  (BJC)  nearest  station  to  the  rapid  terrain  rise  of  the  Front  Range, 
reported  nooderate  snow  between  0055  and  0400  UTC.  Farther  out  in  the  plains, 
Denver  (DEN)  had  periods  of  moderate  snow  between  2219  -  2302  UTC  19  January, 
0030  -  0138  UTC  20  January,  and  0332-0452  UTC. 
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The  period  of  moderate  snow  at  Denver  between  2219  and  2302  appears  related 
to  a  phenomena  other  than  low  level  upslope  flow.  Reflectivity  CAPPIs  at  5  km  in 
Figs.  3.4a-b  show  an  upper  level  band  crossing  over  Denver  during  one  of  the  periods 
of  moderate  snow.  There  was  a  corresponding  region  of  higher  reflectivities  at  lower 
levels  as  well.  The  band  also  had  cores  of  higher  reflectivities  similar  to  cores  described 
by  Carbone  and  Bdine  (1975)  in  their  examination  of  cellular  regions  of  snow 
goieration.  They  found  upper  level  reflectivity  maxima  were  associated  with  vertical 
velocities  of  ±  1.5  ms*^  and  deduced  that  a  trail  of  falling  ice  crystals  associated  with  the 
cores  played  a  role  in  snow  crystal  growth  at  lower  levels,  similar  to  the  "seeder-feeder" 
mechanism  discussed  by  Rutledge  and  Hobbs  (1983). 

Bands  with  embedded  reflectivity  cotes  were  a  common  feature  in  this  case.  The 
video  loop  of  MHR  1.6  elevation  angle  scans  showed  well  deflned  bands  and  cores  in 
the  upper  levels.  In  Figs.  3.4a-d,  a  band  and  associated  cores  identifled  by  the  yellow 
and  orange  tints  representing  higher  reflectivities  illustrate  this  movement.  The  west- 
southwest  to  east-northeast  movement  of  the  cotes  at  speeds  of  approximately  18  ms*^ 
mirrored  the  CLASS  repotted  winds  near  5  km  which  ranged  between  17  and  21  ms'^ 
with  the  stronger  winds  reported  by  the  southern  stations.  Examination  of  the  video 
loop  showed  tiiat  bands  propagated  to  the  southeast  and  at  a  slower  speed  than  the 
cores.  It  appeared  that  the  band  motion  reflected  the  synoptic  scale  motion  associated 
with  the  upper  level  trough  of  low  pressure  which  was  moving  to  the  southeast  about  14 
tns*^  as  estimated  by  the  movement  of  the  500  mb  trough  axis. 

The  0044  UTC  MHR  3.5  and  5  km  (MSL)  reflectivity  CAPPI  in  Figs.  3.5a-b 
show  die  leading  edge  of  another  band  south  of  the  radar.  Vertical  cross  sections  of 
reflectivity  and  radial  velocity  are  shown  in  Figs  3.6a-b.  Note  the  region  of  almost  zero 
velocity  near  18  km  in  the  horizontal  and  5.5  km  in  the  vertical.  A  case  could  be  made 
for  clockwise  rotation  (looking  upwind  towards  the  west  or  counterclockwise  looking 
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downwind  towards  the  east)  around  an  axis  near  20  km  south  of  the  radar.  Hence 
upward  vertical  motion  to  the  south  of  the  axis  and  downward  (or  weaker  upward) 
motion  to  the  north  would  be  expected.  The  reflectivity  cross  section  supports  this 
observation  with  stronger  echoes  to  the  south  in  the  upward  vertical  motion  region  and 
weaker  echoes  in  the  north. 

Another  type  of  banded  or  wavelike  feature  is  shown  in  Fig  3.7.  Here  a  highly 
banded  feature  had  a  series  of  higher  reflecdvity  cores  spaced  at  regular  intervals  of 
approximately  12  km  within  the  band.  Figs.  3.8  a-b  show  cross  sections  of  reflectivity 
and  radial  velocity  taken  through  the  axis  indicated  by  the  black  line  in  Fig.  3.7.  The 
periodic  and  wave-like  appearance  of  the  echoes  suggest  a  different  formation 
mechanism  than  that  which  caused  the  east-northeastward  propagating  cores  in  the 
bands  described  above.  The  periodic  cores  observed  in  this  band  showed  little  motion 
relative  to  the  band.  Satellite  imagery  (not  shown)  indicated  the  band  extended  to  the 
Colorado-Utah  border  which  is  beyond  the  range  of  the  radar,  thus  the  band  appears  to 
end  scmie  70  km  west  of  Pike's  Peak.  This  is  fortuitous  since  it  permits  display  of  the 
terrain  upwind  of  the  band  which  is  dominated  by  the  abrupt  rise  of  the  Continental 
Divide,  illustrated  by  the  gray  shaded  region  in  Fig.  3.7  and  west  of  the  lower  and  flatter 
terrain  of  South  Park.  The  12  km  spacing  between  reflectivity  cores,  lack  of  band 
parallel  core  nnotion  to  the  east-northeast,  and  upwind  terrain  features  are  characteristics 
found  with  trapped  lee  waves  (Duiran  1986).  The  wave  like  structure  in  the  band  was  a 
persistent  feature  and  is  seen  in  Figs  3.4a-d.  Chapter  4  presents  more  discussion  on  the 
possibility  of  mountain  waves  and  their  role  in  precipitation  bands. 

3.3  Extended  Velocity-Azimuth  Display  (EVAD)  Analysis 

Radar  data  can  provide  insight  into  the  kinematic  stmcture  of  mesoscale  features 
beyond  that  obtained  from  simple  examination  of  horizontal  or  vertical  cross  sections. 
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In  general,  dual-Doppler  techniques  provide  an  effective  means  of  determining 
properties  such  as  divergence,  vorticity,  and  horizontal  and  vertical  velocities  on  the 
order  of  a  kilometer  (Srivastava  et  al.  1986)  and  is  limited  by  radar  beamwidth  (which 
increases  with  range),  ambiguities  associated  with  the  intersection  angle  between  the  two 
radar  beams,  and  gridpoint  spacing  (Davies- Jones  1979,  Mohr  et  al.  1986).  Estimates 
of  the  above  properties,  particularly  the  vertical  velocity,  over  a  larger  scale  represents 
another  problem.  Averages  of  the  grid  point  values  over  the  dual-Doppler  domain  could 
be  generated  as  one  way  to  look  at  larger  scale  vertical  motions.  However,  this  method 
is  subject  to  errors  introduced  by  the  interpolation  scheme  as  well  as  errors  in  the 
measurement  of  divergence.  Divergence  errors  may  introduce  errors  into  the  vertical 
motion  calculation  that  are  even  larger  than  the  actual  vertical  velocities  under 
consideration  for  areas  of  widespread  precipitation. 

Early  radar  meteorology  studies  showed  the  potential  capability  of  Doppler  radar 
to  evaluate  the  horizontal  kinematic  structure  of  widespread  precipitation  echoes. 
Lheimitte  and  Atlas  (1961)  laid  the  foundation  for  the  method  used  here  with  the 
velocity-azimuth  display  (VAD)  technique.  Browning  and  Wexler  (1968)  refined  their 
idea  so  titat  VAD  gave  reliable  estimates  of  various  horizontal  kinematic  properties  over 
the  radar  volume  depth.  These  properties  included  divergence,  wind  speed  and 
directimi,  deformation  and  axis  of  dilation,  and  are  the  results  of  calculation  using 
coefficients  obtained  from  a  Fourier  series  analysis  of  the  radial  wind  Held.  Browning 
and  Wexler  (1968)  discussed  several  limitations  to  their  technique  including 
inhomc^eneities  in  precipitation  fall  speed  and  vertical  wind  shear  effects.  They  suggest 
using  elevation  angles  less  than  27  (9)  degrees  in  snow  (rain)  to  avoid  large  errors  in 
VAD  divo’gaice  calculations  introduced  by  the  hydrometer  terminal  fallspeed,  which 
gains  computatitmal  significance  with  increasing  elevation  angles.  The  VAD  technique 
has  been  used  to  evaluate  horizontal  divergence  and  vertical  velocity  with  success  by 
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several  winter  storm  researchers  including  Lilly  (1981),  Shields  et  a/.  (1991),  Hobbs  et 
al.  (1980),  and  Houze  et  al.  (1981). 

The  extended  velocity-azimuth  display  (EVAD)  was  designed  by  Srivastava  et 
al.  (1986)  to  minimize  the  sources  of  error  identified  by  Browning  and  Wexler  (1968). 
The  VAD  technique  is  extended  by  incorporating  the  hydrometeor  terminal  fall  speeds, 
avoided  in  VAD,  into  the  divergence  calculation.  Also,  EVAD  makes  use  of  higher 
elevation  angle  sweeps  than  used  in  VAD  to  obtain  better  estimates  of  terminal  fall 
speeds.  The  increased  accuracy  of  EVAD  divergence  calculations  permit  nxire  accurate 
estimates  of  the  mean  vertical  velocity  incoiporating  all  the  data  in  a  user  specified 
cylindrical  volume  around  the  radar.  This  technique  was  originally  designed  to  study 
the  stratiform  regions  associated  with  mesoscale  convective  systems  (Srivastava  et  al. 
1986).  Matejka  and  Srivastava  (1991)  provide  suggestions  on  data  collection  techniques 
in  which  they  recommend  at  least  20  sweeps  with  a  maximum  elevation  angle  of  at  least 
50  degrees.  In  this  study  the  maximum  elevation  angles  were  22.5  and  17  degrees  and 
the  number  of  scans  in  each  volume  were  24  and  21  for  CSU-CHILL  and  MHR 
respectively.  Since  the  depth  of  the  storm  was  less  than  6  km,  and  horizontal  velocities 
were  fairly  uniform,  it  was  felt  that  the  data  loss  from  lack  of  high  elevation  angles 
would  not  be  too  significant.  Results  from  the  EVAD  technique  appear  reasonable  and 
agree  with  values  for  similar  types  of  storms  over  northeast  Colorado  using  different 
methods  (Reinking  and  Boatman  1986,  Wesley  1991).  However,  the  data  provided  by 
increasing  the  elevation  angle  might  provide  more  accurate  estimates  of  terminal  fall 
speeds  which  would  lead  to  better  divergence  results,  primarily  at  high  elevation  angles 
where  terminal  fall  speeds  gain  significance  in  the  calculation.  An  experiment  to 
compare  true  EVAD  scans  with  the  scanning  strategy  used  during  this  case  would  help 
resolve  ambiguities  over  the  accuracy  of  the  results  presented  below.  EVAD  analyses 
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were  computed  approximately  every  half  hour  and  examined  for  continuity  and 
agreement  between  both  radars. 

The  EVAD  program  is  a  two  step  process,  VAD  analysis  is  run  first  to  obtain 
coefficients  used  in  the  radial  wind  equation  computed  over  a  series  of  rings  along  a 
constant  elevation  angle  scan  from  the  radar  to  the  maximum  cylinder  radius.  The 
coefficients  obtained  from  all  sweeps  constitute  a  data  set  in  which  values  are  identified 
by  range  and  altitude  above  the  ground.  In  EVAD,  a  finite  depth  is  selected  (usually  the 
range  gate  spacing)  over  which  coefficients  havdng  altitudes  within  the  specified  layer 
are  used  to  calculate  the  divergence,  and  terminal  fall  speed.  This  process  is  repeated  for 
each  layer  until  the  top  (bottom,  if  using  top  down  integration)  is  reached.  Once  the 
divCTgence  profile  was  obtained,  the  venical  velocity  was  computed  using  the  anelastic 
version  of  the  continuity  equation; 

Vji'V  +  d  (pw)ldz  =  0  (1) 

Sensitivity  tests  were  conducted  to  examine  the  results  of  varying  the  cylinder 
radius  over  which  the  EVAD  was  calculated.  Within  VAD  and  EVAD,  calculations  of 
wind  coefficients  and  divergence/terminal  fall  speed  can  be  performed  over  user  defined 
radii.  For  the  sensitivity  tests  conducted  in  this  research,  the  largest  radius  used  was  40 
km  in  order  to  avoid  interference  from  mountain  returns,  or  deleted  data  gates  associated 
with  the  mountains.  After  several  iterations  of  various  VAD  and  EVAD  radii 
combinations,  it  was  found  that  the  EVAD  results  were  fairly  insensitive  to  changes  in 
the  VAD  radius  as  long  as  the  VAD  radius  remained  larger  than  the  EVAD  radius  and 
was  greater  than  20  km.  Results  were  most  sensitive  to  the  choice  of  the  EVAD  radius 
particulaily  when  it  was  less  than  20  km.  Increasing  sensitivity  with  decreasing  ring 
size  was  probably  a  result  of  using  a  maximum  elevation  angle  of  22.5  degrees.  The 
sensitivity  study  permitted  optimization  of  the  VAD  ponion  of  the  program  since  more 
con^uter  processing  time  was  needed  for  larger  VAD  radii.  A  VAD  cylinder  diameter 
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of  25  km  was  chosen  as  the  best  compromise  between  processing  time  and  results.  A 
radius  of  24  km  was  chosen  for  the  EVAD  portion  of  the  program  to  take  advantage  of 
as  many  rings  as  possible. 


Results  of  the  temporal  continuity  comparison  were  encouraging,  particularly 
with  the  CSU-CHILL  data.  Figs.  3.9a-d  and  Figs.  3.10a-d  show  EVAD  vertical 
profiles  at  the  same  half  hour  intervals  beginning  2202  UTC  19  January  for  CSU- 
CHILL  and  MHR  respectively.  The  vertical  velocity  profiles  for  the  CSU-CHILL 
EVADs  were  consistent  over  time  and  indicated  a  roughly  parabolic  shape  to  the  profile 
with  a  maxima  near  2800.m  There  is  rather  poor  temporal  agreement  between  MHR 
generated  profiles  at  first  glance.  Upon  closer  inspection,  several  features  persist 
through  the  series  including  peak  upward  motion  near  3000  m.  Profiles  from  MHR  at 
later  times  showed  better  continuity  between  successive  plots  than  those  illustrated  in 
3.10a-d. 

In  both  the  CSU-CHDLL  and  MHR  profiles,  the  straight  line  at  the  top  from  0  to 
the  next  data  point  down  indicates  a  weakness  in  the  procedure.  It  was  assumed  that  the 
vertical  velocity  was  zero  at  the  radar  cloud  top.  This  is  probably  not  the  case  given  the 
insensitivity  of  10  cm  radar  to  small  ice  particles  suspected  at  the  cloud  top.  The  actual 
cloud  top  probably  extends  beyond  the  radar  cloud  top  by  perhaps  as  much  as  1km. 
Soundings  indicate  that  the  tropopause  was  near  7500  m  (MSL)  which  represents  a 
difference  of  at  least  1  km  to  the  highest  echo  tops.  Matejka  and  Srivastava  (1991) 
suggest  that  a  value  other  than  zero  for  the  top  boundary  condition  may  give  better 
results. 

Inter-radar  comparison  shows  agreement  in  the  occurrence  of  a  low  level 
maxima  between  2500  and  3000  m  and  an  poorly  defined  minima  near  4000  m.  Beyond 
that  there  is  not  much  agreement  in  magnitude  or  profile  shape.  Vertical  velocities 
derived  from  the  EVAD  technique  using  CSU-CHILL  appeared  to  have  vertical 
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velocities  approximately  8-10  cms'^  greater  than  results  obtained  from  MHR  data  which 
may  stem  from  two  possibilities: 


1.  CSU-CHILL  had  more  scans  and  a  higher  maximum  elevation  angle  than 
MHR.  This  would  lead  to  more  data  points  for  EVAD  calculations. 

2.  CSU-CHILL  and  MHR  experienced  different  meteorological  conditions. 

Despite  the  differences,  EVAD  analyses  showed  upward  vertical  motion  through  the 
entire  depth  of  the  cloud  for  nearly  ail  profiles.  The  peak  magnitudes  (18  to  25  cms'l) 
from  this  study  compare  favorably  with  observations  and  modeling  studies  of  Front 
Range  snow  storm  vertical  velocities  obtained  by  other  researchers  and  are  shown  in 
Table  1. 


Table  1:  Summary  of  Vertical  Velocities  for  Front  Range  Upslope  Storms. 


Author 

0*orM* 

Vertical 

Velocity 

(cms*l) 

Remarks 

■inryiPiinM— 

'6 

5 

Derived  from  VAD  analysis 

O 

■■ 

Eastern  Wyoming  Anticyclonic 
winter  storm 

BESnHBMIi 

O 

mm 

Shallow  anticyclonic  system, 
terrain  induced 

O 

■■ 

Associated  with  strong  frontal 
passage  at  BAOt 

O 

6.5 

Cyclonic  storm 

Raddatz  and  Khandekar 
(1979) 

M 

15 

Modeling  study  of  upslope 
precip  over  western  Canadian 
Plains 

Abbs  and  Pielke  (1986) 

lO 

Model  results  at  1.35  km 
(AGL) 

■■ 

mm 

Model  study  of  cold  air 
damming 

*  O  =  Observational  study.  M  =  Modeling  study, 
t  BAO  -  Boulder  Atmospheric  Observatory;  3(X)  m  instrumented  tower. 


3.4  Dual'Doppler  Observations 
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One  of  the  earliest  (and  few)  dual-Doppler  studies  of  winter  storms  was  that 
performed  by  Bohne  (1979).  He  examined  shallow  snow  bands  (<  2.5  km  in  depth) 
that  were  associated  with  a  deep,  vertically  stacked  low  pressure  area  over  central 
Illinois.  The  bands  did  not  have  the  well  defined  thermal  wind  orientation  found  in  the 
bands  studied  here.  However,  as  with  this  case,  the  bands  propagated  in  a  manner  that 
suggested  a  connection  to  the  synoptic  flow,  while  precipitation  cores  moved  in  the  band 
with  the  mean  wind  vector  through  the  layer.  Bohne  (1979)  also  used  dual-Doppler 
techniques  to  evaluate  the  vertical  motion  distribution  and  found  a  maximum  magnitude 
of  1.25  ms*^  though  the  mean  magnitude  was  typically  20-40  cms'l.  VAD  analysis 
was  not  performed  on  this  case.  Heymsfield  (1979)  conducted  a  dual-Doppler  study  of 
precipitation  bands  associated  with  a  warm  front  over  northern  Illinois.  Three  major 
bands  were  spaced  about  1(X)  km  apart  and  associated  with  organized  vertical 
circuladons  in  a  layer  2.5  km  deep  just  over  the  warm  frontal  boundary.  Synoptically, 
both  cases  were  different  from  each  other  as  well  as  the  case  examined  here.  Finally, 
Kessinger  and  Lee  (1991)  describe  the  results  of  real-time  operations  of  a  dual-Doppler 
netwoilc  during  WISP-90  strictly  fiom  the  standpoint  of  computer  software  evaluation. 

Dual-Doppler  analysis  in  this  study  was  accomplished  by  using  three  programs 
to  edit,  transfonn  data  from  spherical  to  Cartesian  coordinates,  and  synthesize  data  from 
two  radars  into  a  single  data  set.  As  described  in  section  3. 1 ,  data  were  edited  with 
RDSS.  For  the  Cartesian  transformation,  the  Sorted  Position  Radar  INTerpolation 
(SPRINT)  (Mohr  et  al.  1986)  program  was  used.  The  program  requires  specification  of 
the  X,  y,  and  z  ranges  of  the  Cartesian  coordinate  axis  system  as  part  of  the  program 
input  file.  The  coordinate  ranges  were  selected  on  the  basis  of  data  availability,  crossing 
angle  of  the  radar  beams,  and  the  spatial  resolution.  The  last  two  criteria  are  geometrical 
factors  that  determine  the  accuracy  of  the  dual-Doppler  derived  horizontal  wind  field 
(Davies- Jones  1979).  For  this  case  a  30  degree  crossing  angle  and  1.25  km  spatial 
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resolution  were  used  to  define  the  prime  coverage  area  having  acceptable  horizontal 
velocity  errors.  Usually,  there  are  two  dual-Doppler  lobes,  one  on  each  side  of  the  base 
line  (see  Fig  1.1).  Missing  data  from  MHR  created  large  data  gaps  in  the  northeast 
sector  of  the  volume  precluding  use  of  the  eastern  dual-Doppler  lobe.  The  box  shown  in 
Fig.  1.1  shows  the  region  selected  for  the  horizontal  transform.  Grid  spacing  was  1  km 
in  both  the  x  and  y  directions.  Vertical  spacing  was  500  meters  from  2  to  7  km  MSL 
(0.5  to  5.5  km  AGL).  Figure  3.13  shows  a  close-up  of  the  dual-Doppler  region,  radar 
locations,  the  geometric  discriminators,  and  terrain  contours. 

Once  the  (Cartesian  transform  of  the  radial  velocity  and  reflectivity  fields  from 
both  radars  was  complete,  a  synthesis  was  performed  using  the  Custom  Editing  and 
Display  of  Reduced  Information  in  Cartesian  Space  (CEDRIC)  program  (Mohr  et  al. 
1986)  to  obtain  u  (east-west)  and  v  (north-south)  winds,  convergence  Helds,  and 
vertical  velocity.  Vertical  velocities  were  obtained  by  integration  of  the  anelastic  form  of 
the  continuity  equation,  assuming  w  =  0  at  the  top  of  each  echo  column.  The  two  radar 
solution  obtained  in  C!EDRIC  necessitated  computing  hydrometer  fallspeeds  from  a 
reflectivity-fallspeed  relationship  derived  by  Atlas  et  al.  (1973). 

Unfortunately,  the  first  attempt  at  a  dual-Doppler  analysis  indicated  there  was  a 
problem  in  the  chain  of  programs  leading  to  the  dual-Doppler  synthesis.  There  were 
unacceptable  g2q)S  in  the  vertical  velocity  field  across  a  horizontal  plane,  furthermore, 
gaps  increased  with  height.  After  many  sensitivity  tests  in  which  input  and  output 
parameters  for  RDSS,  SPRINT,  and  CEDRIC  were  altered  in  a  variety  of  ways  for  both 
radars  to  find  maximum  coverage,  the  conclusion  was  reached  that  a  problem  existed 
with  the  SPRINT  interpolation  scheme.  The  problem  was  particularly  noticeable  with 
C^U-CHILL  data  and  may  be  a  result  of  fewer  data  points  from  running  the  radar  at 
slightly  reduced  power  due  to  transmitter  arcing.  An  independent  test  run  of  CSU- 
CHDLL  data  by  CSU-CHILL  personnel  using  other  transformation  software  developed 
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at  the  University  of  Chicago,  verified  SPRINT  was  not  inteipolating  all  the  RDSS 
universal  format  spherical  coordinate  data  to  Canesian  coordinates.  It  is  unclear  where 
the  problem  with  SPRINT  lies  and  time  constraints  prohibited  an  examination  of  the 
software  or  seeking  an  alternative  program. 

Despite  the  problems  encountered  in  generating  vertical  velocities  frcrni 
synthesized  fields,  enough  data  was  extracted  to  prepare  horizontal  wind  vector  fields  at 
most  levels  for  the  2328  UTC  19  January  volume  scan.  Figs  3.14a-f  show  the  MHR 
CAPPI  reflectivity  patterns  and  horizontal  winds  obtained  from  the  dual-Doppler 
analysis.  An  arrow  with  the  length  of  one  gridpoint  spacing  represents  a  speed  of  15 
ms*^  Comparison  between  dual-Doppler  derived  winds  and  the  Denver  wind  profiler 
(Rg.  2.7b)  shows  good  agreement.  Low-level  flow  was  north-northeasterly  backing 
with  height,  then  increasing  in  speed  with  little  change  in  direction  above  4  km.  Note 
the  weak  winds  at  3  km  (Fig.  3.14c),  this  corresponds  to  the  layer  of  maximum  shear 
between  the  northeasterly  surface  winds  associated  with  the  arctic  airmass  and  the  west- 
southwesterly  winds  aloft.  Tables  2  and  3  list  the  mean,  standard  deviation,  minimum, 
and  maximum  of  the  u  and  v  wind  components  computed  for  each  level  between  2  and  5 
kmMSL. 

Finally,  a  bulk  Richardson  number  profile  was  computed  using  the  mean  wind 
values  fen-  the  dual-Doppler  re^on  and  thermodynamic  data  from  the  (XX)0  UTC  20 
January  Berthoud  sounding.  The  bulk  Richardson  number  was  calculated  using  Stull’s 
(1988)  equation; 

R  g  AOyAz 

®~ev[(AU)2  +  (AV)2] 

where : 

g=s  gravity 

AOy  s  the  change  in  the  virtual  potential  temperature  over  a  layer 

Az  =  Layer  depth 
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6v  =  Mean  virtual  potential  temperature  in  the  layer 

AU  =  Change  in  the  u  (east- west)  wind  component  over  the  layer 

AV  =  Change  in  the  v  (north-south)  wind  component  over  the  layer 


The  results  of  the  calculation  are  summarized  in  Table  4.  There  appears  to  be  good 
agreement  between  the  dual-Doppler  derived  Richardson  number  and  the  sounding 
derived  values  which  are  discussed  in  Chapter  4. 


Table  2:  Dual-Doppler  u  Wind  Component  Statistics. 


■  Helgh't  ■ 
(km) 

Mean  u 
(ms'l) 

u  Std 
Dev 
(ms'i) 

4.7 

rm— 

wmm 

^.5 

msmm 

miiam 

1  3.5 

mMm 

-$.3 

msam 

wwm 

msism 

i£2Qiai 

-4.8 

mmm 

mmimm 

-5.3 

m^im 

-6.1 

mssm 

HESSI 

wssxm 

3.i 

wsim 

1  6.5 

maim 

4.7 
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Table  3:  Dual-Doppler  v  Wind  Component  Statistics. 


Height 

(km) 

Mean  V 
(ms*l) 

V  Std 
Dev 
(ms-1) 

V  Min 
(ms-1) 

vMax 

(ms-1) 

2.0 

-7.20 

1.38 

-11.5 

-0.9 

2.5 

-11.2 

“Tin 

575 

"^4:55 

-l7.7 

“SIT? 

3.5 

-1.45 

“05 

9.1 

4.0 

~r55 

2.03 

-7.7 

”30 

4.5 

~Tl5 

”05 

“TTO 

1.01 

"or 

-14.2 

5.5 

5.5 

1.51 

”353 

“TO 

1.S3 

-ili 

“Tin 

6.5 

3.32 

“430 

”353 

7.7 

Table  4:  Dual-Doppler  Richardson  Number  Calculadon 


Center  of  Layer 

Richardson  Number 

i.i5 

l5.55 

03 

535 

5.55 

053 

533 

1.57 

433 

1534 

4.75 

55.01 

333 

534 

5.75 

41.83 

533 

TO! 

3^  Summary 

Radar  data  collection  and  prcx^ssing  techniques  were  discussed  in  this  chapter  as 
wdl  as  single-  and  dual-Dq)pler  analyses.  Results  of  EVAD  calculations  of  vertical 
vekxdty  were  also  presented.  Single-Doppler  analysis  showed  the  banded  nature  of 
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echoes  associated  with  this  storm  including  some  wave-like  bands  over  the  southern 
portions  of  the  radar  coverage  area.  Higher  reflectivity  cores  were  embedded  within  the 
bands.  The  cores  propagated  with  the  upper  level  winds.  Vertical  velocity  profiles 
derived  firom  the  E  VAD  technique  compared  favorably  with  values  obtained  by  other 
researchers.  Problons  with  the  interpolation  program  precluded  calculations  of  the 
vertical  velocity  distribution  using  dual-Doppler  analysis.  However,  good  agreement 
was  found  between  horizontal  winds  computed  from  dual-Doppler  analysis  and  the 
DEN  profiler.  The  next  ch^ter  will  address  possible  mechanisms  for  snow  band 
finmation. 
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a)  MHR  REHECnVITY  (dB^  2249  UTC  -  25  KM 
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Hgure  3.1:  2249 ITTC  19  January  Mile  High  Radar  (MHR)  constant  altitude  plan 
mtion  indicator  (CAPPD  display  for  a)  23  km  MSL  and  b)  3.5  km  MSL.  Color  bar  at 
boom  indicaiBS  ranges  of  reflectivity  in  dBZ  for  each  shading,  east-west  and  north- 
soath  distances  are  in  km.  Gaps  in  coverage  from  the  southwest  to  ncnthwest  are  due  to 
ddetion  of  mountain  returns. 


Figure  3.  Ic-d:  Same  as  (a)  except  for  c)  4.5  km  MSL  and  d)  5.5  km  MSL. 
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Hsure  3^2:  2249  UTC  19  January  MUe  High  Radar  CAPPI  display  of  radial  velocity 
for  ^  2.5jan  MSL  and  b)  3.5  km  MSL.  Color  bar  at  bottom  indicates  ranges  of  radial 
vdocity  in  ms'^  for  each  shading,  east-west  and  north-south  distances  are  m  km.  Dashed 
ccmtours  indicate  negative  velocities  (towards  the  radar),  solid  contours  represent 
po^tive  velocities  (away  from  the  radar).  Coverage  gaps  to  west  are  due  to  mountains. 
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a)  MHR  RETlECrnvnY  (dB^  2348  UTC  -  25  KM 
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b)  miR  ranirnvnY  (dE2^  0044  urc  -  25  km 


-40 


-00 


-00  -40  -ao  0  ao  40  00 

Dhtann  EMt-ltat  (nO 

0  8  10  a*  00 

Figure  3.3:  2.5  km  MSL  Mile  High  Radar  CAPPI  reflectivity  display  of  reflectivity  for 
a)  2346  UTC  19  January  and  b)  0044  UTC  20  January.  Color  bar  at  bottom  indicates 
ranges  of  reflectivity  in  dBZ,  cast- west  and  north-south  distances  are  in  km.  Gaps  in 
coverage  from  the  southwest  to  northwest  are  due  to  deletion  of  mountain  returns. 
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MHR  REFLBCTIVrrY  (dBZ)  2211  UTC  -  5  KM 
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Figure  3.4:  5  km  constant  altitude  plan  position  indicator  CAPPI  reflectivity  display 
for  a)  221 1  UTC  and  b)  2230  UTC  19  January  Color  bar  at  bottom  indicates  ranges 
of  reflectivity  in  dBZ,  east- west  and  north-south  distances  are  in  km.  Gaps  in  coverage 
from  the  southwest  to  northwest  are  due  to  deletion  of  mountain  returns. 
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a)  MHR  RETlECnVITY  (dBZ)  0044  UTC  -  aS  KM 
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Figure  3.5:  0044  UTC  20  January  CAPPI  reflectivity  display  for  a)  3.5  km  MSL  and 
b)  5.0  km  MSL.  Color  bar  at  bottom  indicates  ranges  of  reflectivity  in  dB2^  east-west 
and  north-south  distances  are  in  km.  Gaps  in  coverage  to  west  are  due  to  deletion  of 
mountain  returns.  North-south  line  shows  cross  section  transect  used  for  Figs.  3.6a-b. 
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a)  railBCrriVlTY  (dBZ)  CROSS  SECTION  X  =  0 
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Figure  3.6:  (X)44  UTC  20  January  Mile  High  Radar  vertical  cross  section  of  a) 
reflectivity  and  b)  radial  velocity.  Color  bar  at  bottom  indicates  ranges  of  reflectivity 
in  dBZ.  South  is  to  the  left.  West  is  looking  into  the  figure.  Dashed  lines  in  Fig  3.6b 
indicate  regions  of  negative  (towards  the  radar)  velocity. 
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MHR  REFLECTIVITY  (dBZ)  2249  UTC  - 
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Rgmc  3.7:  2249  UTC  19  January  Mile  High  Radar  CAPPI  display  of  reflectivity  for  4.5 
km  MSL.  Red  lines  are  tmain  contours  in  intervals  of 200  m  to  2000  meters  then  every 
500  m.  Dark  gray  indicates  terrain  between  3000  and  3500  m.  Li^t  gray  sha^ng 

rqnesents  toiain  above  3500  m.  APA  indicates  locaticm  of  Centennial  airport  Color  bar  at 
bottom  indicates  ranges  ci  reflectivity  in  dBZ  for  each  shading.  Black  line  represents  axis 
used  in  Hgs.  3.8a'b. 
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Hgiffc  3.8:  2249  UTC  19  January  Mile  High  Radar  vertical  cross  section  of  a)  reflectivity 
md  b)  radial  vdochy.  Color  bar  at  bottom  indicates  ranges  of  reflectivity  in  dBZ  for  each 
riiflding,  Horiztmtal  distance  is  measured  from  the  left  end  of  the  black  line  in  Fig.  3.7. 
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Figure  3.9:  CSU-CHILL  EVAD  vertical  velocity  profile  for  a)  2202  UTC,  b)  223 1  UTC 
c)  2259  UTC,  and  d)  2328  UTC  19  January. 
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Figure  3.1 1:  Dual-Doppler  synthesis  region.  CSU-CHILL  and  MHR  are  indicated  by 
CHL  and  MHR.  Laige  area  within  arc  represents  regions  where  crossing  angle  is 
greater  dian  30  Agrees.  Area  within  curv^  triangle  represents  region  where  the  spatial 
resoluticMi  is  less  than  1.25  km.  Terrain  contours  are  shown  for  every  200  m  to  2000  m 
and  dark  and  light  gray  shading  highlight  terrain  between  2000-3000  and  3000-4000m 
respectively. 
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Dual  Dc^plo*  Analysis 


Hgure  3.12  (a):  2.0  km  (MSL)  2328  UTC  19  January  Dual-Doppler  derived  wirid  vectors 
and  MHR  reflectivity  field.  Vector  length  of  cme  gridpoint  spacing  represents  15  ms'^ 
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Figure  3.12  (b):  Same  as  (a)  except  for  2.5  km. 
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Hgure  3.12  (c):  Same  as  (a)  except  for  3.0  km. 
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Dual  Dc^pler  Analj^  2236  UTC  -  3^ 


Figure  3.12  (d):  Same  as  (a)  except  for  3.5  km. 
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Figure  3.12  (e):  Same  as  (a)  except  for  4.0  km. 
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Figure  3.12  (g):  Same  as  (a)  except  for  5.0  km. 


108 


Chapter  4 


POSSIBLE  SNOW  BAND  FORMATION  MECHANISMS 

As  discussed  in  Chapter  1,  mesoscale  precipitation  bands  have  been  observed  in 
diverse  geographical  locations  and  under  a  wide  range  of  synoptic  situations.  Various 
dynamic  mechanisms  have  been  proposed  to  explain  the  formation  of  these  bands 
including:  craiditicMial  and  craivective  instabilities,  conditional  symmetric  instability 
(CSI),  internal  and  ducted  gravity  waves,  Kelvin-Helmholtz  waves,  gravity  currents, 
and  upper  level  jet  streaks.  However,  terrain  influences  caused  by  the  Rocky  Mountains 
add  mechanisms  not  sufficiently  addressed  in  other  major  winter  storm  studies.  Wesley 
(1991)  discusses  the  role  of  a  barrier  jet  in  causing  heavy  snow  along  the  Front  Range  in 
cases  of  strong  cold  air  damming.  Precipitation  bands  over  eastern  Colorado  may  also 
be  die  result  of  mountain  wave  activity  (Wesley  and  Pielke  1990).  In  this  chapter, 
various  mechanisms  are  examined  for  their  potential  in  forming  the  snow  bands 
observed  during  this  study. 

4.1  Conditional  Instability 

Instability  as  a  result  of  con'  we  processes  in  the  vicinity  of  frontal  ztHies  can 

cause  [recantation  bands  (Snook  1992),  similar  to  mechanisms  responsible  for  squall 
line  activity  (Hane,  1986).  Emanuel  (1983b)  menticms  that  the  susceptibility  of  the 
atmoqdiere  to  moist  adiabatic  vertical  motions  from  conditiraial  instability  governs  the 
structure  of  bands  associated  with  CSI.  Conditional  instability  was  assessed  using 
Stull's  (1988)  method  to  evaluate  its  role  in  creating  bands  or  influencing  CSI 
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circulations.  StuU  (1988)  states  that  air  can  be  considered  conditionally  unstable  when 
the  dry  static  stability  gradient  is  positive  (ds/dz  >  0)  and  the  saturation  static  stability  is 
negative  (9Ses/3z  <  0).  Where: 

S  =  C^T  +  gz;  Dry  static  energy  or  Montgomery  stream  function 
Ses  =  S  +  Lyrgat:  Saturation  static  energy 

Cp  =  1004.67  m2  s'2  K*l;  Specific  heat  of  dry  air  at  constant  pressure 
Lv  *  2.5  X  10®  J  kg*f ;  Latent  heat  of  vaporization 
rsat  -  Saturation  mixing  ration  for  water 

Vertical  profiles  of  s  and  Ses  change  were  prepared  fex  0000  UTC  20  January 
soundings  and  examined  for  regions  of  conditional  instability.  Figs.  4.  la-b  show 
vertical  profiles  dry  and  saturated  static  energy  change  (ds/dz  and  dSes/dz  respectively) 
at  Akron  and  Elbert,  CO.  Conditional  instability  criteria  were  met  where  the  dashed  line 
(saturatkm  static  energy  change)  is  negative  and  the  solid  line  (dry  static  stability  change) 
is  positive.  The  Akron  data  shows  conditional  instability  was  limited  to  a  layer  between 
5000  and  6500  MSL  while  the  Elbert  profile  indicates  a  more  unstable  environment  with 
lower  over  all  values  compared  to  Akron.  Akron  was  located  in  a  region  of  diminished 
»iow  band  activity  con^ared  to  the  pronounced  bands  in  the  vicinity  of  Elbert  The 
Elbert  profile  also  shows  shallow,  low  level  layers  of  conditional  instability  within  1000 
m  of  the  sui^ce  (below  2500  m  MSL)  and  a  distinct  inclination  towards  conditional 
instability  above  4000  m  MSL  where  well  defined  bands  were  observed  in  the  radar 
data.  Low  level  convergence  (and  upward  votical  motion)  in  the  upslope  flow  over  the 
northern  portiems  of  the  Palmo’  Divide  may  have  played  a  role  in  the  release  of 
conditional  instability  and  subsequent  CSI  snow  band  formation  (Emanuel  1983a,  b).  It 
is  unlikely  that  the  shallow  layers  of  conditional  instability  near  the  surface  alone  are 
reqxHisible  for  snow  band  development 

4.2  Conditional  Symmetric  Instability 
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Conditimal  symmetric  instability  (CSI)  and  its  role  in  the  formation  of 
precipitation  bands  has  been  studied  by  many  of  the  winter  storm  researchers  cited  in 
this  thesis.  Bennetts  and  Hoskins  (1979)  proposed  precipitation  bands  formed  as  a 
result  of  CSI  in  a  layer.  Emanuel  (1983a,  b)  simplified  the  process  by  introduction  of  a 
Lagrangian  parcel  method  of  CSI  evaluation  which  has  since  been  embraced  by  most 
researchers.  Briefly,  CSI  theory  states  that  slantwise  convection  is  the  result  of 
combining  convective  and  inertial  instabilities  that  arise  from  an  unstable  distribution  of 
body  forces  acting  on  parcels  within  a  fluid.  Cbnvective  instabilities  cause  vertical 
motKHi  since  the  body  force  is  gravity,  whereas  inertial  instability  causes  radial  motion 
on  a  horizontal  plane  where  centrifugal  body  forces  dominate  (Emanuel  1983a). 
Symmetric  instability,  the  sum  of  these  two  instabilities,  produces  slanted  motion  with 
respect  to  tiie  horizontal  and  vertical  directions. 

In  the  case  of ‘^purely  meridional,  steady,  moist  baroclinic  shear  flow  in  thermal 
wind  balance”,  Emanuel  (1983  a,b)  defines  a  quantity  Mg,  the  geostrophic  pseudo- 
angular  ooomentum  by: 


MgSVg  +  fx, 

(3) 


tHdiere  is  the  geostrophic  flow  normal  to  an  x-z  cross  section  of  Mg  surfaces,  /is  the 
Coricdis  parameter,  x  is  the  distance  along  the  jc  axis  of  the  plane,  increasing  firom  left  to 
tight  If  the  pressure  gradient  in  the  y  direction  is  zero,  then  the  y  component  of  the 
momentum  equation  can  be  written  as: 


dvg 

dt 


(4) 
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indicating  that  Mg  is  conserved  following  an  air  parcel.  From  Emanuel  (1983b),  the 
horizontal  (x)  and  vertical  (z)  components  of  the  momentum  equation  are: 


^  =/(v  -Vg)^f{Mp-Mg),  (5) 

-  Bvg),  (6) 

Ovo 

where  By  is  the  virtual  potential  temperature,  subscript  p  refers  to  the  parcel  M  or  6^, 
siinilaily  g  refers  to  geostn^hically  balanced  M  or  of  the  free  atmosphere,  o  refers 
to  a  reference  6^.  Acceleration  of  an  air  parcel  in  the  positive  (negative)  x  directicHi  takes 
place  if  (Mp  -  Mg)  in  Eq.  3  is  positive  (negative),  implying  the  parcel  nx>mentum  is 
greater  (less )  than  the  environmental  momentum.  Vertical  acceleration  of  an  air  parcel 
tq)  (down)  will  occur  if  (6^  -  6vg)  is  positive  (negative)  indicating  a  parcel  virtual 
potential  tenqrerature  greater  (less)  than  the  environment  Snook  (1992)  points  out  the 
requirement  for  slantwise  convection  is  met  when  {Mp  <  Mg)  and  (Byp  >  Byg).  From 
die  inacdcal  stanc^nt  of  evaluating  CSI  on  vertical  cross  sections,  diis  implies  Mg 
must  decrease  in  the  positive  x  direction  along  a  ctmstant  Be  line  (Sanders  and  Bosan 
198Sa,  Wolfsberg  et  al.  1986).  There  are  several  quick  qualitative  assessments  that  will 
he4>  in  evaluating  the  potential  for  CSI.  From  Snodt  (1992)  and  Seltzo*  et  al.  (1986); 

1)  Winds  increase  with  height  indicating  a  baroclinic  atmosphere. 

2)  A  unidirectional  wind  profile  (though  not  always  necessary). 

3)  A  neaiiy  saturated,  well-mixed  vertical  thermodynamic  profile  with  neutral  to 
near  moist  adiabatic  Itqise  rate. 

4)  Observations  of  multiple  bands.  Firontogenetic  forcing  can  «cplain  single 
bands,  Init  not  multiple  bands. 

5)  Bands  orknted  parallel  to  the  thermal  wind. 
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6)  A  region  in  the  atmosphere  where  the  Richardson  number  is  less  than  one. 

Using  the  qualitative  indicators  above,  data  were  checked  fen*  CSI  potential 
before  preparatitxi  of  a  CSI  cross  section.  Examination  of  wind  profiles  (see  for 
example  Figs.  2.7a-b  and  2.10a-c)  showed  that  winds  had  an  almost  unidirectional 
profile,  and  increased  with  height  (above  the  arctic  airmass),  which  can  be  taken  as 
evidence  of  the  required  atmospheric  baroclinicity.  A  nearly  saturated  atmosphere  with 
kq)se  rates  tqqnoaching  moist  adiabatic  was  indicated  by  the  0000  UTC  20  January 
soundinp.  Multiple  bands  were  observed  by  radar,  particularly  over  the  southern 
portitxis  of  die  project  area.  Precipitation  bands  visible  in  the  radar  data  nearly  paralleled 
the  thermal  wind  (compare  Figs.  3.4a-d  with  Fig.  2.4). 

The  requirement  for  a  Richardson  number  less  than  one  occurred  at  several 
levds  (Fip.  4.2a-b).  The  Richardson  numbers  obtained  from  the  dual-Doppler 
syndiesis  (Table  4)  ccxnpare  favorably  with  numbers  obtained  from  soundings.  The 
radar  derived  values  have  the  advantage  of  being  calculated  with  an  average  horizontal 
wind.  The  disadvantage  of  this  technique  comes  from  the  use  of  relatively  thick  shear 
layers.  Stull  (1988)  cautions  that  Richardson  number  calculations  frxxn  soundings  (even 
more  so  w4ien  using  the  coarse  vertical  resoluticm  in  the  radar  data  analysis  technique) 
are  prone  to  inaocurades  due  in  part  to  ^proximation  of  the  vertical  partial  derivatives 
by  finite  differences.  It  is  possible  that  actual  Richardsem  numbers  were  lower  than 
those  calculated,  particulaily  if  there  were  large  gradients  over  a  shaUow  depth.  If  this 
hypothesis  is  true,  tiien  the  dual-Doppler  drived  values  might  have  been  lower  than 
those  in  Table  4.  In  summary  however,  both  sounding  and  dual-Doppler  Richardson 
number  calculatkms  met  CSI  criteria. 

Once  it  was  established  diat  tire  qualitative  criteria  for  CSI  were  met,  a  cross 
section  of  Afg  and  0g  contours  was  prepared.  The  vertical  cross  section  was  taken  along 
die  line  shown  on  Hg.  2.4.  This  is  approximately  orthogonal  to  the  geostrophic  wind  at 
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nxKt  levels.  Mg  ivas  calculated  using  Eq.  3,  then  plotted  and  analyzed.  The  analysis 
was  combined  with  the  equivalent  potential  temperature  analysis  and  is  shown  in  Fig. 
4.3.  Solid  lines  represent  in  *K  and  dashed  lines  represent  Mg  in  ms‘^  Following 
CSI  convention,  regions  of  possible  CSI  are  indicated  by  the  hatched  region  where  Mg 
decreases  along  a  constant  de  line.  Two  regions  of  possible  CSI  are  apparent,  one  near 
the  surface  and  the  other  in  the  middle  troposphere.  The  region  near  the  surface  should 
be  igntned  as  the  6*e  analysis  is  suspect  at  this  level  and  is  probably  an  artifact  of  the 
analysis  routine,  since  surface  pressures  are  close  to  850  mb  along  the  cross  section. 
Howevo-,  at  upper  levels,  CSI  could  have  played  a  role  in  forming  bands.  There  is 
good  agreement  witii  radar  observations  sintre  they  show  most  bands  were  well  defined 
above  4  km  (~615  mb).  However,  despite  the  favorable  indications  of  CSI,  other 
mechanisms  still  must  be  ccmsidered  for  their  role  in  snow  band  development 

4,3  Ducted  Gravity  Waves,  Internal  Gravity  Waves,  and  Kelvin* 
Helmholtz  Instabilities 

Ltndzen  and  Tung  (1976)  proposed  ducted  gravity  waves  as  another  cause  of 
precipitation  bands.  Originally  directed  towards  explaining  bands  associated  with 
convective  activity,  others  have  applied  the  theory  to  winter  storms  (Parsons  and  Hobbs 
1983,  Wolfsberg  era/.  1986).  Parsons  and  Hd}bs  (1983)  examined  internal  gravity 
waves  and  ducts  as  a  system  in  relation  to  g^erating  and  maintaining  warm  sector  and 
wide  ctdd-fttMital  rain  bands.  In  the  case  of  warm  sector  rainbands,  they  found  some 
agreement  between  theory  and  obsovation,  however  low  level  potential  instability  and 
lack  of  a  stable  duct  forced  them  to  conclude  that  this  rainband  type  was  associated  witii 
odier  mechanisms.  They  assessed  gravity  wave  potential  as  being  even  lower  fn*  the 
wide  cttid-fiontal  rain  band  due  to  a  critical  level  in  the  duct,  and  slow  phase  speed  of 
the  band  cooqiared  to  predicted  propagation  speeds  of  waves  in  a  low  Richardson 
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number  regime.  Wolfsberg  et  al.  ( 1986)  noted  the  vertical  atmospheric  profile  in  the 
case  diey  studied  was  favraable  for  a  gravity  wave  duct  but  they  could  find  no  evidence 
of  a  surface  pressure  perturbation  and  thus  concluded  bands  were  the  result  of  other 
mechanisms. 

In  the  case  studied  here,  the  possibility  of  internal  and  ducted  gravity  waves  was 
considered  and  rejected  (except  internal  waves  associated  with  mountain  waves 
discussed  in  section  4.4)  as  a  cause  of  major  snow  bands  fcH-  the  following  reasons: 

1.  Bands  were  not  orthogonal  to  the  mean  flow  as  required  by  ducted  gravity 
wave  theory,  but  instead  were  parallel  to  the  mean  flow. 

2.  Examination  of  surface  pressure  traces  showed  no  evidence  of  perturbations 
which  Lindzen  and  Tung  (1976)  use  as  a  measure  of  “response”  to  ducted  waves. 

3.  Location  of  the  bands  with  respect  to  synoptic  fronts  did  not  favor 
pn^gation  of  waves  fhxn  frontal  disturbances.  Bands  occurred  well  behind  the 
surface  cold  front,  unlike  the  cases  examined  by  Parsons  and  Hobbs  (1983). 

The  presence  of  Kelvin-Helmholtz  (K-H)  instabilities  as  a  cause  for  major  band 
formation  was  discounted  since  again,  the  OTthogcmal  flow  requirement  (bands  were 
parallel  to  the  flow)  was  not  met  (Stull  1988).  However,  within  the  bands,  flow 
condidiMis  could  be  favorable  for  K-H  waves  transverse  to  the  major  bands. 

Richardson  number  proBles  at  Elbert  and  Denver  (Hgs.  4.2a-b)  show  shallow,  well 
defined  regions  of  low  Richardson  numbers  (less  than  one)  just  above  SOOO  m  in  both 
proBles  which  were  similar  to  profiles  from  other  soundings  (not  shown).  Stull  (1988) 
states  that  in  the  presence  of  laminar  flow  the  onset  of  K-H  instability  waves  requires  a 
lUchardson  number  <  0.2S.  However,  if  the  atmosphere  is  already  turbulent, 
Richardson  numbers  <  1  indicate  turbulent  conditions  will  persist  and  K-H  instabilities 
can  occur  until  the  Richardson  number  increases  to  greater  than  one.  Stull  (1988) 
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remaiks  tfiat  K-H  waves  can  have  horizontal  scales  up  to  hundreds  of  kilometers  but  are 
limited  in  vertical  extent  to  tens  to  hundred  of  meters.  Richardson  number  profiles 
indicate  favcnable  conditions  for  K-H  wave  formation  in  this  storm,  however  there  was 
no  means  to  confirm  or  deny  their  existence. 

4.4  Gravity  Currents 

Hobbs  and  Persson  (1982)  used  the  concept  of  gravity  currents  as  a  possible 
explanation  for  narrow  cold-frontal  bands  observed  during  the  CYCLES  program. 
Convective  instabilities  arise  from  bulges  and  clefts  associated  with  a  denser  fluid 
advancing  into  a  region  of  lower  densides.  Bulges  at  the  top  of  the  denser  fluid  come 
fitxn  internal  circulations  in  the  fluid.  Oefts  are  the  result  of  frictional  effects  cm  the 
denser  fluid  by  the  lower  boundary.  A  nose  of  denser  fluid  protrudes  over  a  shallow 
layer  of  less  dense  fluid  near  the  lower  boundary.  This  causes  ctxivective  instabilities 
and  bands  near  the  frontal  boundary  (Parsons  and  Hobbs  1983).  The  cold  front  in  the 
case  examined  here  had  passed  though  the  region  several  hours  before  weU  wganized 
banded  features  were  observed,  consequently  a  gravity  current  mechanism  for  snow 
band  fonnation  is  not  ccxisidered  likely  since  this  phenomena  is  expected  in  the  vicinity 
of  thefrxHit 

4.5  Upper  Level  Jet  Streaks 

Dunn  (1988)  showed  the  importance  of  ctxisidering  upward  vertical  moticms 
associated  with  jet  streaks  when  evaluating  the  potential  for  snow  band  formation  over 
northeast  Colorado.  Uccellini  and  Kocin  (1987)  examined  the  role  of  jet  streaks  in 
producing  heavy  snow  over  the  Eastern  Seaboard.  Upward  vertical  motions  associated 
widi  jet  streaks  under  no  tempoature  advection  conditions  typically  occur  in  favored 
quadrants  which  are  the  left  exit  and  right  entrance  regims  when  viewed  looking  down 
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wind  (Keyser  and  Shapiro  1986).  Fig.  4.4a  illustrates  the  no  temperature  advection 
example  ctf  a  jet  streak  where  isotherms  parallel  die  axis  of  the  jet  streak.  Mid- 
trc^sphoic  vertical  nnotion  maxima  (illustrated  by  +  and  in  pressure  coordinates)  are 
located  symmetrically  about  the  jet  streak  axis.  One  would  expect  to  find  regions  of 
enhanced  precipitation  in  those  areas  under  the  negative  sign  (i.e.  left  exit  and  right 
entrance).  Keyser  and  Shapiro  (1986)  also  examined  the  effects  of  temperature 
advection  on  the  location  of  vertical  motion  centers  with  respect  to  the  jet  streak.  Hg. 
4.4b  is  dieir  example  of  a  warm  air  advection  case  and  is  similar  to  conditions  observed 
over  northeast  Colorado  (compare  with  Fig.  2.5c).  Note  the  northward  shift  of  the 
upward  vertical  motion  maximum  in  the  entry  region  as  a  result  of  warm  air  advection. 

Examination  of  the  NMC  300  mb  analysis  (Fig.  2.Sc),  Platteville  wind  proHler 
data  (Fig.  7.a),  and  time  cross  sections  (Hgs  2.13a-d)  shows  jet  stream  level  winds 
decreased  (i.e.  the  core  of  the  jet  streak  has  moved  east)  and  warm  air  advection  at  jet 
stream  level.  Furthermore,  closer  inspection  of  cross  sections  (Figs  2.12a  and  2.13a-d) 
and  300  mb  CLASS  staticm  plots  (Figs.  2.6a-d)  show  the  axis  of  the  jet  was  near  the 
southern  stations.  If  the  region  of  vertical  velocity  was  shifted  northward  ^m 
tempoature  advection,  the  southern  portions  of  the  project  area  would  experience  greater 
upward  motion  and  precipitation.  Radar  data  suppons  this  conclusion  since,  in  general, 
bands  with  higher  radar  tops  were  observed  to  the  south.  Also,  there  was  a  region  of 
higher  snowfall  over  the  Palmer  Divide  (see  Fig.  2.1)  that  ctmesponded  with  radar 
obsmed  reflectivity  bands.  Finally,  satellite  imagery  (not  shown)  showed  the  major 
band  over  the  Palmer  Divide  extended  west  beyond  the  radar  coverage  area  to  the  axis  of 
the  300  mb  trough  near  the  Utah-Colorado  border. 

4.6  Terrain-Induced  Mechanisms 
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During  arctic  outbreaks  the  rising  terrain  of  the  Front  Range  acts  as  a  banier  to 
the  cold  air  and  effectively  blocks  the  flow  from  further  westward  movement  (Wesley 
1991).  Furthermore,  in  some  cases  the  forced  lifting  of  warmer,  moist  air  over  the  cold 
stable  air  below  can  result  in  damming.  Wesley  and  Pielke  (1990)  showed  examples  of 
blocking  and  damming  and  discussed  how  these  phenomena  can  enhance  precipitation 
over  die  eastern  plains  of  Ck)lorado  adjacent  to  the  Front  Range.  In  the  case  of 
damming,  snow  bands  form  along  low  level  convergence  lines  that  develop  near  the 
Ftcmt  Range  in  response  to  cyclonic  turning  of  the  easterly  flow  due  to  its  inability  to 
ascend  higher  terrain.  Snow  bands  that  develop  during  blocking  conditions  are  often 
pmiallel  to  the  thermal  wind.  There  is  little  correlation  between  bands  and  surface 
omvergence  zones  except  in  extreme  cases  of  blocking  when  development  of  a  low  level 
banier  jet  fnms  parallel  to  the  Front  Range.  This  can  lead  to  a  band  of  heavier  snow 
just  east  of  the  Front  Range  (Wesley  1991).  The  existence  of  damming  and  location  of 
convergence  zones  can  easily  be  determined  by  checking  isotherm  and  streamline 
analyses.  Surface  temperatures  will  decrease  from  east  to  west  and  the  streamline 
analysis  should  indicate  rapid  cyclonic  turning  in  the  easterly  flow  near  the  Front  Range. 
Isotherm  analyses  (Hg.  2. 16  for  example)  for  the  case  under  study  here  showed  no  well 
defined  cold  pool  near  the  mountains,  nor  did  stream  line  analyses  (Fig.  2.17  for 
example)  show  any  easterly  flow,  let  alone  any  well  defined  convergence  zones 
associated  with  cyclonic  turning.  Therefore,  we  conclude  that  cold  air  damming  was  not 
present  in  this  case.  While  blocking  conditions  were  present  as  illustrated  by  the 
inability  of  die  cold  air  to  ascend  over  the  mountains  and  west  southwesterly  flow  aloft, 
there  was  no  evidence  of  a  low  level  barrier  jet  in  (XASS  or  wind  profiler  data.  Snow 
bands  aligned  themselves  with  the  thermal  wind,  which  was  similar  to  the  blocking  case 
discussed  by  Wesley  and  Pielke  (1990). 
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One  of  the  most  striking  features  in  Figs.  3.1a-d  and  3.4a-d  is  the  band  south  of 
the  radar  having  rather  evenly  spaced  reflectivity  cores  with  an  approximate  interval  of 
12  km.  Fig.  3.7  and  Figs.  3.8a-b  provide  a  more  detailed  look  at  the  horizontal  and 
vertical  structure  of  the  band.  Observations  showed  there  was  little  movement  of  the 
cores  within  the  band,  unlike  the  single  precipitation  cores  discussed  in  Chapter  3.  The 
cores  in  this  band  appeared  to  develop  and  dissipate  with  little  down  wind  movement,  in 
a  band  relative  sense.  The  periodic  nature  of  the  reflectivity  pattern  suggests  wave 
features  were  embedded  in  the  the  larger  scale  band.  Durran  (1986)  states  trapped 
mountain  waves  can  form  in  the  lee  of  a  mountain  barrier  if  certain  conditions  are  met. 
Fig.  4.5  shows  a  model  of  the  air  flow  associated  with  trapped  mountain  waves.  These 
waves  are  the  result  of  upward  propagating  waves,  initiated  by  the  barrier,  being 
reflected  downward  by  an  upper  boundary  (a  stable  layer)  which  in  turn  become  upward 
propagating  waves  again  on  striking  the  ground.  In  the  ideal  case  there  will  be  no  loss 
in  wave  amplitude  since  there  is  no  energy  exchanged  during  the  reflections,  thus  the 
process  will  (theoretically)  create  an  infinite  number  of  waves.  The  visible  downstream 
waves,  called  trapped  waves,  are  the  superposition  of  the  upward  and  downward 
propagating  waves  and  typically  have  wavelengths  between  5  and  25  km  (Durran  1986). 

Calculation  of  the  Scorer  parameter  (Scorer  1949)  represents  a  relatively  simple 
means  oi  evaluating  the  potential  for  trapped  mountain  waves.  The  Scorer  parameter, 
can  be  calculated  from: 


r2 

~U2  U  dz2 


(7) 


whoe  is  the  Brunt-VflisHlS  frequency,  and  U  is  the  wind  speed  at  the  center  of  the 
layer  under  evaluation.  In  practice,  Durran  (1986)  recommends  excluding  the  second 
tenn  finm  die  equation,  as  it  is  difficult  to  evaluate  from  the  coarse  vertical  resolution  of 
rawinsonde  observations  and  its  contribution  to  the  result  is  small  compared  to  the  first 
term.  If  die  Scorer  parameter  decreases  significandy  with  height,  then  the  potential  for 
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trapped  waves  exists.  Queney  et  al.  (1960)  identified  other  qualitative  requirements  for 
strong  lee  wave  development: 

1.  The  mountain  barrier  must  have  a  steep  lee  slope.  The  strongest  waves  occur 
when  there  is  a  gradual  windward  slope  and  steep  lee  slope. 

2.  The  wind  must  have  a  direction  within  30*  of  normal  to  the  ridge  line  through 
a  deep  layer.  Ridge  top  wind  speed  must  be  at  least  7-15  ms-l  (depending  on  the 
character  of  the  ridge),  and  wind  speeds  must  increase  with  height  above  the  ridge. 

3.  An  inversion  or  stable  layer  must  be  present  at  the  ridge  top  level  and 
upstream  of  the  barrier  with  weaker  stability  aloft 

Examination  of  Grand  Junction,  CX)  sounding  (Fig.  4.6)  shows  the  top  of  the 
low  level  inversion  at  672  mb  or  3310  m  with  decreasing  stability  aloft  The  sounding 
also  shows  winds  were  within  the  30*  of  perpendicular  to  nnountain  tops  if  the  axis  of 
the  Continental  Divide  is  assumed  to  be  oriented  north-south  over  central  Colorado. 
Furthermore,  winds  were  greater  than  13  ms*l  at  mountain  top  level  and  increased  with 
height  Scorer  parameter  calculations  were  performed  for  all  CLASS  stations,  Denver, 
and  Grand  Junction.  Grand  Junction,  Denver,  Elbert,  and  Berthoud  (Figs.  4.7a-d) 
Scorer  parameter  profiles  had  many  similarities  including  a  maximum  near  3.5  km 
(MSL)  (except  Grand  Junction),  and  a  rapid  decrease  to  a  minimum  near  6  km  (MSL). 
Hgure  3.7  shows  the  waves  were  organized  downstream  of  the  Continental  Divide, 
which  has  a  steeper  lee  slope  compared  to  its  windward  slope. 

Scorer  (1949)  devel(q)ed  an  equation  to  calculate  the  wave  number  of  tipped 
waves  in  a  two  layer  atmosphere.  According  to  Durian  (1986)  use  of  his  equation  can 
jvovide  useful  estimates  of  the  wavelength.  From  Scorer,  the  wave  number,  k,  of  the 
tnqjped  wave  can  be  obtained  from  the  following  expression: 

(L;2-  k  2)*/2  COt[(L/ -  k  2)*'2//]  +  (it  2  _  £^2)1/2  ^  Q  (g) 
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Where: 


Lj  =  Scorer  parameter  in  the  lower  layer 
L2  =  Scorer  parameter  in  the  upper  layer 
H  =  Depth  of  the  lower  layer 
k  -  Wave  number  of  the  trapped  wave. 

Values  for  and  were  obtained  from  the  lowest  two  layers  in  the  Elbert  Scorer 

parameter  profile  (Fig.  4.7c)  where  there  was  the  smallest  change  with  height.  Thus  the 
mean  values  of  and  £2^  were  3.162  X  10'^  and  3.162x  10'^  m'^  respectively 

corresponding  approximately  to  +0.5  and  -0.5  on  the  x-axis  of  figure  4.7c.  H  was 
estimated  to  be  1000  m.  By  numerical  methods,  k  was  found  to  be  5.845  x  10^  m'^ 
which  equates  to  a  wavelength  (A.)  of  approximately  10.7  km  since  X  =  hUk.  This 
numerically  derived  value  of  the  wavelength  compares  favorably  with  the  12  km  value 
obtained  fiom  direct  observation  which  lends  suppon  to  the  conclusion  that  trapped 
waves  were  present. 

Determination  of  trapped  wave  amplitudes  is  a  complex  problem  with  no  simple 
solutitm  as  it  depends  on  the  interaction  of  several  factors:  the  scale  and  shape  of  the 
upwind  banio’,  vertical  wind  stratification,  and  vertical  stability  (mainly  temperature) 
stratification.  Coiby  and  Wallington  (1956)  developed  a  relationship  that  anempts  to 
explain  the  vertical  displacement  of  a  streamline  that  incorporates  the  factors  listed  above 
fiom  which  one  might  infer  characteristics  about  the  wave  amplitude.  Their  equation  is: 

C,  =  -2nihbe-kb)  {UJU,)  WzjcidWojJdkY^  sinfct  (9) 

A  B  C 


Where: 

^2  =  Votical  displacement  of  a  streamline 
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h  =  Height  of  barrier 
b  =  Half  width  parameter  of  the  barrier 
k  =  Wave  number  of  the  trapped  wave 
U  =  Horizontal  wind 

Y  =  Stability  parameter,  satisfies  the  condition: 

+  (L^  +  k^)Y =0  (10) 

L  =  Scorer  parameter 

X  =  Horizontal  distance 

Subscripts  o,  z,  and  k  represent  with  respect  to  the  surface,  some  height 
z,  and  wave  number  respectively. 

A,  B,  and  C  represent  the  terms  associated  with  moditications  of  the  streamline 
displacement  as  a  result  of  barrier  characteristics,  wind  stratification,  and  stability 
stratification  respectively. 

De^te  the  difficulties  encountered  in  computing  the  amplitude  of  the  waves,  Corby  and 
Wallington  offo’  scmie  qualitative  observations  and  comments.  Rom  (9),  the  wave 
an^litude  dqrends  in  a  complex  manner  on  the  topography  and  the  vertical  distribution 
of  atmospheric  prcqrerties.  However,  in  examining  the  elements  of  the  equation,  two 
statements  can  be  made  about  the  amplitude; 

1.  The  maximum  amplitude  will  occur  when  the  barrier  width  is  approximately 
equal  to  the  airstream  wavelength  (i.e.  -kb  =  0). 

2.  Large  amplitude  waves  occur  when  a  shallow  layer  of  great  stability  rather 
than  a  deep  layer  of  slight  stability  is  present  (i.e.  term  C  is  large). 

Non  linearities  make  solution  of  (9)  extremely  difficult  so  that  at  best  only 
genoel  statements  can  be  made  about  the  wave  amplitude.  Unfortunately,  no  means  of 
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d^ermining  the  amplitude  of  the  waves  from  observations  (hence  inferences  of  wave 
induced  vertical  velocities)  was  available  during  the  experiment 


The  enhanced  reflectivity  cores  in  the  waves  (see  Figs.  3.7  and  3.8a)  probably 
represent  regions  where  upward  ascent  in  the  flow  modified  the  ice  crystal  number 
density  and/or  radar  back  scatter  cross  section.  Similarly,  weaker  reflectivities  between 
the  cores  probably  correspond  to  regions  of  descent  Radar  theory  (see  Battan,  1973) 
states  that  the  reflectivity  factor  Z  (or  equivalent  reflectivity  factor  2^  in  this  case)  for  a 
given  volume  is  directly  proportional  to  the  summation  of  scatterer  number  density,  n, 
of  diameter,  D,  times  the  diameter  raised  to  the  sixth  power  or: 

Z=  SnjD^mm^m"^  (11) 

i 


Also, 


dBZ=101ogZ  (12) 

Examination  of  Fig  4.8a  shows  that  the  maximum  difference  between  cores  and  weak 
echo  regions  was  approximately  25  dbZ.  Thus  AZ  would  be  lO^-^  or  316  mm^  m'^. 
Since  temperatures  w^e  below  freezing  through  the  atmosphere,  assume  that  nearly  all 
the  radar  energy  was  scattered  by  ice  crystals.  If  the  ice  crystal  diameters  were  all  the 
same,  die  number  concentration  would  have  to  increase  (decrease)  cm  the  order  of  300  to 
account  for  the  increase  (decrease)  in  returned  power  associated  with  the  cotes  (weak 
echo  regions). 

A  more  likely  cause  for  the  observed  structure  in  the  radar  echo  is  a  change  in  the 
ice  crystal  diameter.  If  the  number  concentration  stayed  constant  and  all  ice  crystals 
were  of  the  same  initial  diameter,  then  an  increase  (decrease)  in  the  diameter  by  a  factor 
of  2.6  would  explain  the  increase  (decrease)  in  returned  power.  If  the  level  of  maximum 
ascent  rate  cxteurred  in  a  regicm  where  the  temperature  regime  was  between  -12  and 
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>16*C  ice  crystals  would  experience  maximum  growth  rates  according  to  Auer  and 
White  (1982).  Consider  an  ice  crystal  in  water  saturated  air  at  a  temperature  between  -5 
and  -20  *C  traveling  alcaig  the  rising  half  of  a  sinusoidal  wave  that  is  9  km  long  at  a 
speed  of  25  ms'^  Under  these  conditions  the  ice  crystal  will  grow  by  deposition  during 
an  ascent  that  would  last  for  360  seconds  (6  minutes),  neglecting  the  increase  in  its 
terminal  velocity  as  its  mass  increases.  The  ice  crystal  should  increase  in  diameter  by  at 
least  a  factor  of  2.6  during  the  dme  of  ascent  As  the  air  descends  during  the  other  half 
of  the  cycle,  adiabatic  warming  and  lower  ice  saturation  may  cause  the  crystal  to 
sublimate,  dius  reducing  the  amount  of  radar  energy  scattered  back  to  the  receiver. 
Growdi  and  decay  of  ice  crystals,  rather  than  large  change  in  number  density,  seems  a 
much  more  reasonable  explanation  for  the  observed  echo  pattern. 

Since  it  is  unlikely  that  a  unifcsm  size  distribution  existed,  the  actual  cause  for 
die  periodic  structure  in  the  radar  echoes  was  probably  due  to  a  complex  interaction 
between  deposition  and  sublimation,  aggregation,  riming  as  crystals  travel  through 
tegkms  of  high  liquid  water  content,  precipitation  of  ice  crystals,  and  changes  in  the 
number  density.  Detailed  in-situ  measurements  of  microphysical  and  kinematic 
properties  throughout  the  volume  of  the  wave  clouds  would  have  provided  much  greater 
insight  into  the  processes  involved  in  creating  the  observed  periodic  structure. 
UnfortunatBly,  profiler  data  fnxn  the  Palmer  Divide  region  were  not  available.  Aircraft 
data  was  limited,  however,  a  pass  over  the  Palmer  Divide  by  the  University  of 
Wyoming  King  Air  after  0040  UTC  at  3500  m  MSL  showed  the  traiperature  to  be 
between  -10  and  -12  'C.  Crystal  habits  included  a  mixnne  of  graupel  and  dendrites, 
including  aggregates  and  rimed  crystals.  Dendrites  seemed  to  predominate  indicating 
vapor  depodtion  was  the  primary  cause  of  crystal  growth,  at  least  in  the  upper  levels. 
Riming  probably  aided  in  growth,  particularly  in  the  lower  levels  whote  higher  liquid 
water  ctmtents  were  likely.  Radiometer  data  from  Elbert  indicated  periodic  depletions  of 
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the  vertical  integrated  liquid  water  content  between  2200  UTC 19  January  and  0000  20 
January.  Values  ranged  between  0  and  0.44  mm,  while  the  integrated  water  vapor 
varied  between  0.S0  and  0.67  cm.  Greater  snowfall  (or  at  least  higher  liquid  water 
equivalent)  was  likely  during  periods  of  high  liquid  water  content  since  rimed  crystals 
with  greater  terminal  velocities  and  aggregates  would  have  been  present  These  periodic 
fluctuations  in  the  radiometer  data  were  probably  associated  with  the  movement  of  other 
features  in  the  bands  and  not  the  periodic  reflectivity  cores  since  they  showed  little 
movement  with  respect  to  the  band. 

The  question  of  how  non-trapped  wave  precipitation  cores  were  formed  still 
remains  largely  unanswered.  Wesley  and  Helke  (1990)  suggested  that  droplets  forming 
in  nucleation  zones  upwind  of  rtK>untain  peaks  first  reported  by  Rauber  et  al.  (1986), 
mi^t  undergo  freezing  and  sufficient  lifting  to  clear  the  barrier  in  regions  of  mountain 
wave  vertical  ascent  Alternatively,  ice  nucleation  from  deposition  or  freezing  may  take 
place  in  the  wave  above  the  mountain  barrier.  Rogers  and  Vali  (1987)  found  that  clouds 
in  contact  with  a  snow  surface  experienced  an  increase  in  ice  nucleation  perhaps  due  to 
ice  crystals  carried  aloft  by  the  airflow.  Regardless  of  the  primary  mechanism, 
seccMidary  ice  formaticm  processes  such  as  splintering,  fragmentation  of  crystals,  and 
ctMitact  freezing  as  discussed  by  Hobbs  and  Rangno  (1985)  would  multiply  the  number 
of  ice  nuclei. 

The  shallow  convective  cell  theory  proposed  by  Carbone  and  Bohne  (1975)  and 
a  similar  seeder-feeder  model  developed  by  Rutledge  and  Hobbs  (1983, 1984) 
represents  anodier  aiqnoach  where  seeder  cells  dqxrsit  ice  crystals  into  feeder  clouds  of 
relatively  higher  liquid  water  content.  The  low  integrated  liquid  wato*  contents 
associated  witii  heavier  snow  may  represent  those  locations  in  a  snow  band  where  liquid 
waio'  in  die  feecto  cloud  has  been  converted  to  ice  crystals  by  seeding  from  uppo*  level 
clouds. 
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Mulvey  (1977)  repented  the  results  of  a  modeling  study  on  the  survival  time  and 
transport  distance  for  ice  crystals  as  functions  of  ice  saturation,  wind  velocity, 
tonperature,  and  ice  crystal  size.  For  example,  an  ice  crystal  of  diameter  1 .7  mm 
initially  fonned  at  -25  'C  and  ice  saturation  of  0.9,  and  experiencing  a  wind  speed 
regime  20  ms'^  had  a  maximum  range  of  347  km.  A  more  realistic  example  was  a 
maximum  distance  of  124  km  experienced  by  an  ice  crystal  of  2.5  mm  diameter  initially 
formed  in  an  ice  saturation  of  0.9  at  -15  *C  and  transported  by  20  ms'l  winds. 

Transport  distance  of  ice  crystals  in  a  subsaturated  environment  is  limited  by  sublimation 
thus  large  crystals  will  have  the  greatest  range.  In  a  super  saturated  environment, 
tranqxxt  distance  is  limited  by  the  fall  speed  of  the  crystal,  therefore  small  crystals  will 
have  the  greatest  range.  Inonporating  the  results  of  Mulvey  *s  study  with  the  Hndings  of 
Rogers  and  Vali  (1987),  under  the  right  saturation  and  temperature  regime,  strong  winds 
at  mountain  level  could  cause  considerable  downwind  transport  of  ice  crystals 
picked  up  from  the  surface.  These  ice  crystals  could  serve  as  nuclei,  or  have  a  role  in  a 
seeder-feeder  mechanism. 

4.7  Discussion  and  Conceptual  Model 

Analysis  of  the  total  snowfall  (Fig.  2.1)  showed  three  distinct  maxima.  What 
were  the  mechanisms  responsible  fex-  the  observed  snowfall  distribution?  There  is  no 
question  tiiat  orographic  lifting  in  conjunction  with  cooling  after  frontal  passage  over 
northeast  Colorado  helped  create  a  nearly  saturated  environment  at  the  surface  and 
saturated  conditions  aloft  Strong  upslq)e  ascent  appears  to  have  occurred  near  Fort 
Collins  and  offers  a  partial  explanation  for  the  maximum  west  of  the  city.  However, 
radar  data  showed  that  an  upper  level  band  crossed  over  the  region  during  the  storm 
which  teads  to  qjeculation  that  upslope  induced  snowfall  was  enhanced  by  a  seeder- 


126 


feedm  mechanism.  Ice  crystals  falling  from  the  band  probably  seeded  the  relatively  high 
liquid  water  content  in  lower  level  clouds. 

The  same  orographic  lifting  and  seeder-feeder  enhancement  couplet  spears  to 
have  occuned  over  the  region  northeast  of  Denver,  though  cnographic  effects  were 
probably  less  inqxsrtant  than  in  the  Fort  Collins  region.  The  seeder-feeder  mechanism 
played  a  mcxie  important  role.  Radar  data  showed  two  distinct  flow  regimes  and 
enhanced  reflectivity  regions.  Within  the  low  level  arctic  airmass  a  well  defined  region 
moved  from  the  north  northeast  with  the  low  level  flow.  The  region  increased  in  areal 
coverage  and  intensified  somewhat  as  it  approached  the  western  edge  of  the  Palmer 
Divide.  Above  the  arctic  airmass  were  well  defined  bands  with  embedded  enhanced 
reflectivity  cores.  The  cores  moved  from  south  southwest  to  north  northeast  Seeding 
of  die  low  level  clouds  by  ice  crystals  may  have  played  a  mrae  important  role  here  than 
over  the  regicm  west  of  Fort  Collins.  As  discussed  in  section  2.3  terrain  induced  vertical 
motkms  were  on  the  order  of  5  cm  s*l,  thus  a  seeder-feeder  mechanism  might  might 
provide  an  explanation  for  the  heavier  snowfall  northeast  of  Denver. 

The  region  of  higher  snowfall  over  the  eastern  portions  of  the  Palmer  Divide 
qjpears  related  to  a  different  mechanism.  Radar  data  and  satellite  imagery  showed  a 
well  defined  ccmtinuous  band  crossing  the  region.  The  band  was  oriented  parallel  to  the 
upper  level  winds  and  was  near  the  northern  edge  of  the  jet  axis.  This  situation  may 
have  been  similar  to  that  examined  by  Dunn  (1988).  In  Dunn's  case,  jet  streak  induced 
vertical  motions  in  a  region  of  conditional  symmetric  instability  were  argued  to  trigger 
slantwise  convection  and  heavier  precipitation.  Radar  data  did  not  show  any  well 
defined  low  kvel  reflectivity  enhancements  (except  the  band)  as  with  the  otiier  two 
regions  described  above. 

Analysis  of  CSI  showed  there  was  potential  for  slantwise  convection  and  thus 
formation  of  bands.  Two  mechanisms  seemed  to  have  played  a  role  in  triggering  the 
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convection,  release  of  conditional  instability,  possibly  from  a  seeder-feeder  couplet,  and 
jet  streak  induced  vertical  motions.  It  would  seem  likely  that  the  conditional  instability 
release  mechanism  might  involve  a  positive  feedback  circuit,  in  which  ice  crystals  falling 
fixxn  a  band  into  the  feeder  zone  release  instability  that  in  turn  leads  to  slantwise 
C(Mivectk>n  which  results  in  more  ice  crystal  precipitation  and  so  on. 

The  process  of  band  formation  is  complex  and  involves  the  interacticMi  of  several 
phoiomena.  However,  a  conceptual  model  of  the  kinematic  structure  is  presented  in 
Fig.  4.8a'b.  Figure  4.8a  shows  a  plan  view  of  a  band  with  an  embedded  transient  core. 
For  die  sake  of  clarity  low  level  clouds  are  not  shown.  The  thin  solid  arrow  parallel  to 
die  band  axis  ^ows  die  band  relative  motion  of  the  core.  The  long  arrows  parallel  to 
the  band  show  the  upper  level  stream  lines.  Short  arrows  pointing  to  the  lower  right 
represent  motion  of  the  band.  Note  that  the  core  moves  parallel  to  the  streamline  and 
bands  move  in  an  off  axis  directirm  that  seems  to  correspond  to  the  movement  of 
synoptic  scale  features.  Dashed  anows  represent  low  level  wind  flow.  The  shear 
region  between  die  low  level  northeasteriies  in  die  arctic  air  mass  and  the  upper  level 
west  southwesterlies  is  shallow,  on  the  order  of  500  to  1000  m. 

Hgure  4.8b  shows  a  cross  section  of  the  band.  Scallt^ied  regions  represent 
clouds.  The  Rocky  Mountains  are  the  higher  terrain  to  die  lefL  The  shear  zone  is  the 
hatched  region  mid  way  up  into  the  cloud.  Winds  are  from  the  north  to  northeast  below 
the  shear  zone,  winds  are  from  the  southwest  to  west  above.  Cores  and  ice  crystals 
falling  into  die  lower  clouds  (the  feeder  zone)  are  shown.  Ice  crystals  in  the  cores 
probably  form  contact  nucleation  as  described  in  Hobbs  and  Rangno  (1985)  then  fall 
into  tire  feetter  zone  where  drey  can  grow  by  deposition,  aggregation,  and  riming,  or 
play  a  rede  in  ice  multiplication.  The  larger  ice  particles  can  then  fall  out  as  snow  (ex* 
graupel  if  heavily  rimed). 

4.8  Summary 
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In  this  chapter  various  mechanisms  were  examined  for  the  role,  if  any,  they 
might  have  played  in  producing  the  upper  level  snow  bands  and  precipitation  cotes 
described  in  the  radar  observations.  Evaluation  of  the  data  suggests  that  conditional 
instability,  jet  streaks,  and  mountain  waves  played  a  role  in  band  formation  and  that  CSI 
analysis  showed  the  atmosphere  was  susceptible  to  slantwise  convection.  Conditional 
instability  probably  played  a  role  in  triggering  bands  over  the  northern  portions  of  the 
project  area,  while  jet  streak  induced  vertical  nnotions  were  nnore  important  to  the  south. 
The  most  striking  examples  of  banding  occurred  over  the  southern  portions  of  the 
project  area  where  it  t^pears  mountain  wave  activity,  in  OMijunction  with  jet  streak 
induced  vertical  motions,  may  have  enhanced  bands. 
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a)  AKR  DHTA  ENHWY  -  0000  UTC  20  JAN 


OHJA  DlBRCy 

b)  HB  Hm  ENERGY  -  0000  UTC  20  JAN 


-e  0  2  4  8  B  10 

ULTA  Bioer  0(**V!9**2) 

4.1  0000  UTC  20  January  dry  and  saturation  static  energy  change  for  a)  Akitm  and  b) 
Berdiond.  Sdid  line  is  diy  static  energy  change,  dashed  line  is  saturation  energy 
change. 
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RICHARDSON  NUMBER 


RICHARDSON  NUMBER 

42  0000  UTC  20  January  Richardson  number  proflle  for  a)  Elbm  and  b)  Denver 
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4.3  Vertical  cross  section  of  Og  and  Mg .  Solid  lines  are  dg  contours  in  *K,  dashed  lines 
are  Mg  contours  in  ms'^  FCL,  DEN,  PD  represent  approximate  positions  of  Fort 
Collins,  Denver,  and  the  Palmer  Divide  respectively. 


a) 


b) 


4.4  Schematic  illustratioii  of  a  jet  streak,  a)  no  temperature  advection  case  and  b)  warm 
air  advection  case.  Thick  solid  lines,  thick  dashed  lines,  thin  solid  lines,  thick  solid 
allows,  and  plus  and  minus  signs  represent  geopotential  height  contours,  isotachs, 
isodierms,  sense  of  cross-flow  ageostrophic  wind  component  at  level  of  maximum 
wind,  and  sense  of  midtropo^hoic  prossure-coordinate  vertical  velocity,  O) 
req)ectively  (after  Keyser  and  Shapiro  1986) 


HEIGHT  (km) 


CROSS-MOUNTAIN  DISTANCE  (km) 


4.5  Model  of  traroed  mountain  waves  in  the  x-z  plane.  Thin  lines  are  streamlines  (after 
Durran  1986). 
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u-f  I  I  r  I - 1  I  I  I 

-‘lO  -30  -20  -JO  0  10  20  30  HO 

TMPC  DHPC 


it*  Junction,  CO  0000  UTC  20  January  rawinsondc  sounding.  Heavy  solid  line 

curve.  Horirontal  wlid  lines 

m  isobars  (mbXwlid  (hagonal  hnes  sloping  from  lower  left  to  uppo-  right  arc 

1 2’tf^  adiabate  m  the  solid  curved  lines  sloping  from  the  lower  right  to 
left,  and  the  moist  adiabats  arc  the  dashed  lines  slewing  from  lower  right  to  upper 

left.  Pennants  on  wind  shafts  represent  25  ms-i,  long  barbs  5  ms-i,  and  short  barbs  2  5 
ms*i. 
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kfio  (aOOBBt  PARAMEIQl  *  IObCO 


faglO  CSOQRER  PARAMETBl  •  lOd^ 


4.7  (a)  0000  UTC  20  Januaiy  Scorer  parameter  profiles  for  a)  Grand  Junction,  b) 
Denver,  c)  Elbert,  and  d)  Beithoud. 
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48 

text  for  cxplanauon. 


a)  plim  view  and  b)  vertical  ptoffl* 
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Chapter  5 


DISCUSSION  AND  CONCLUSIONS 

This  research  used  the  WISP-91  data  set  to  examine  a  winter  storm  that  occurred 
over  northeast  Golcarado  on  19-20  January  1991  with  the  objectives  of  detailing 
observations  of  snow  bands  and  gaining  insight  into  possible  band  formation 
mechanisms.  Radar  data  were  used  to  identify  and  study  the  organization  and  structure 
of  the  bands  as  well  as  provide  quantitative  insight  into  their  kinematic  structure. 

Terrain  features  had  a  signiiScant  impact  on  snow  fall  amounts  over  the  region, 
particttlaily  vk^en  considering  the  low  level  flow  and  its  orientation  to  topographic 
bairiors.  However,  above  the  arctic  air  mass  well  defined  bands  of  higher  radar 
reflectivity  had  similarities  to  precipitatitm  bands  studied  by  other  researchers  with  the 
exceptitm  of  mountain  wave  induced  bands  and  associated  trapped  wave  phenomena. 
Various  causes  for  precipitation  band  formation  published  in  the  literature  were 
evaluated  fw  their  potential  in  causing  the  snow  bands  observed  in  this  case. 

The  WISP-91  data  set  permitted  mesoscale  analysis  of  surface  and  upper  air  data 
obtained  from  59  automated  weather  stations  over  northeast  Colorado,  Wyoming  and 
Nelnaska,  in  addition  to  the  standard  surface  airway  observation  reporting  stations.  The 
data  set  also  included  radar  data  from  CSU  CHILL  and  MHR  radars,  3  hourly  CLASS 
rawinscmde  soundings  from  five  to  six  sites  over  northeast  Colorado,  two  wind 
profilers,  microwave  radiometers,  and  a  special  snow  spotter  network.  Software  was 
writt«i  to  give  existing  meteorological  analysis  and  display  packages  the  capability  of 
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processing  the  volume  of  data  collected  during  the  experiment.  Also,  teirain  data  was 
incorporated  into  analyses,  including  overlaying  radar  data  on  terrain  maps. 

On  the  synoptic  scale,  an  upper  level  trough  of  low  pressure  and  associated 
arctic  cold  front  moved  over  northeast  Colorado.  Post-frontal  northeast  surface  flow 
brought  upslope  conditions  to  the  region  with  up  to  15  cm  (6  in)  of  snow  over  an  18 
hour  period.  The  heaviest  snow  fall  was  organized  into  three  distinct  regions  (see  Fig. 
2.1).  Isotherm  and  streamline  analysis  showed  that  cold  air  damming  as  defined  by 
Wesley  (1991)  was  not  present  in  this  case.  Blocking  of  the  arctic  air  by  the  mountains 
did  occur  but  not  with  sufficient  intensity  to  cause  formation  of  a  low  level  barrier  jet 
StitMig  west-southwesterly  flow  (over  50  ms*^)  aloft  associated  with  a  jet  streak,  and  a 
favorable  thennodynamic  profile,  caused  trapped  mountain  waves  down  wind  of  the 
Qmtinental  Divide.  The  vertical  thermodynamic  profile  also  showed  the  temperature 
was  between  -12  and  -15  *C  near  600  mb  which  is  a  favorable  condition  fw  heavy 
snows. 

Radar  data  showed  mesoscale  snow  bands  were  a  common  feature  of  this  steam. 
There  was  a  preference  for  mwe  intense  band  formation  over  the  southern  portions  of 
the  project  area.  A  video  loop  of  radar  data  that  included  terrain  contours  helped  identify 
two  distinct  flow  features,  a  low  level  regime  that  was  influenced  by  terrain,  and  an 
upper  level  regime  where  bands  moved  in  response  to  fiee  atmosphere  motions.  Vfithin 
the  bands,  higher  reflectivity  cores,  (typically  5-10  dBZ  greater  than  the  surrounding 
band  echo)  were  present.  The  cores  transited  the  band  following  the  upper  level  winds, 
while  the  band  itself  drifted  slowly  to  the  south-southeast.  Radar  observations  of 
trqjped  mountain  waves  illustrated  the  importance  of  considering  mountain  induced 
flow  modifications  in  snowstorms  and  snow  band  formation  over  eastern  Colorado. 

Analysis  of  EVAD  derived  vertical  motion  profiles  from  both  CSU-CHILL  and 
MHR  radar  data,  indicated  similar  values  to  vertical  motion  results  for  northeast 
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Colorado  winter  storms  obtained  by  other  researchers.  The  results  of  the  experiment 
indicated  that  peak  mesoscale  vertical  motions  were  on  the  order  of  0.2  ms'^  between 
2500  and  3000  m  (MSL)  and  were  associated  with  an  ill-deHned  snow  band  and 
associated  region  of  higher  reflectivities  near  the  surface. 

Dual-Doppler  analysis  at  one  particular  time  provided  insight  into  the  horizontal 
flow  structure.  Low  level  northeasterly  winds  backed  with  height  to  a  shallow  shear 
zone  between  3000  and  3500  m  (MSL).  Richardson  number  calculations,  using  mean 
winds  derived  from  the  dual-Doppler  synthesis  region  and  the  Berthoud  Qass  sounding 
thermodynamic  data,  agreed  reasonably  well  with  values  obtained  from  soundings 
alone. 

Of  tiie  band  formation  mechanisms  examined,  conditional  instability  and  jet 
streak  induced  vertical  velocities  could  have  triggered  slantwise  craivection  and  thus 
bands  since  the  atmosphere  showed  potential  for  CSI.  Mountain  waves  also  played  a 
role  in  explaining  the  radar  echo  structure.  Midtropo^heric  vertical  motions  associated 
with  the  entrance  region  of  a  jet  streak  were  shifted  northward  by  warm  air  advectirm 
along  the  jet  and  helped  explain  the  well  defined  band  over  the  Palmer  Divide.  Satellite 
imagery  showed  this  band  extended  west  beyond  the  radar  coverage  area  to  the  axis  of 
the  300  mb  trough  near  the  Utah-Colorado  border.  Other  bands  farther  north  did  not 
appear  to  have  as  strong  a  link  to  jet  streak  vertical  motions,  but  instead  seem  related  to 
ctxiditional  instability  and  CSI  induced  circulations.  Scorer  parameter  calculatitms 
indicated  nearly  stationary  trapped  mountain  waves,  with  wavelengths  of  about  12  km, 
occurred  in  the  band  over  the  Palmer  Divide.  Regions  of  vertical  ascent  associated  with 
the  waves,  secondary  ice  production  methods,  and  transport  of  surface  based  ice 
crystals,  may  have  been  responsible  fcx*  creating  small  ice  particles  which  in  turn  could 
have  been  important  in  a  seeder-feeder  process  (where  clouds  associated  with  low  level 
flow  may  provide  the  feeder  portion  of  the  couplet).  Reflectivity  cores  within  other 
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bands  appeared  to  be  weak  convective  regions,  especially  at  higher  levels,  that  probably 
{ffovided  seed  ice  crystals  in  a  seeder-feeder  process  involving  lower  level  clouds. 

5.1  Suggestions  for  Future  Research 

Several  questicms  could  not  be  completely  answered  in  this  study  and  may  foim 
the  basis  for  further  research.  These  questions  are; 

1.  What  is  the  formation  mechanism  for  the  reflectivity  cores?  Are  they  truly 
associated  with  shallow  convection,  the  result  of  mountain  induced  flows,  or  a 
axnbination  of  both? 

2.  Are  well  defined  bands  necessarily  associated  with  jet  streaks  or  was  this 
case  unique  in  that  respect? 

3.  Are  bands  induced  by  mountain  waves  a  common  feature  and  if  so  are  there 
any  fav(»ed  geographical  regions  for  their  formation? 

4.  What  are  the  detailed  microphysical  processes  associated  with  the  trapped 
wave  phenomena? 

A  radar  based  climatological  study  of  snow  bands  might  help  resolve  the  issue  of 
band  location  with  respect  to  jet  streaks  and  mountain  barriers,  while  aircraft  studies 
would  provide  more  insight  into  the  microphysics  of  bands  associated  with  trapped 
waves.  The  EVAD  technique  looks  promising  as  a  means  of  determining  mesoscale 
vertical  motirms.  A  sensitivity  study  comparing  true  EVAD  volume  scans  with  the 
scanning  strategy  employed  during  this  experiment  would  help  remove  ambiguities 
surrounding  results  obtained  from  low  elevation  angles  used  during  shallow  (compared 
to  thunderstorms)  winter  storms. 


141 


REFERENCES 


Agee,  E.  M.  and  S.  R.  Gilbert,  1989:  An  aircraft  investigation  of  mesoscale  convection 
over  Lake  Michigan  during  the  10  January  1984  cold  air  outbreak.  J.  Atmos. 
ScL,  46,  1877-1897. 

Atlas,  D.,  R.  C.  Srivastava,  and  R.  S.  Sekhon,  1973:  Doppler  characteristics  of 

precipitation  at  vertical  incidence.  Rev.  of  Geophys.  and  Space  Phys.  2, 1-35. 

Auer,  A.  R  Jr.  and  J.M  White,  1982:  The  combined  role  of  kinematics, 

thermodynamics,  and  cloud  physics  associated  with  heavy  snowfall  episodes.  J. 
Meteor.  Soc.  Japan,  60,  500-507. 

Battan,  L.  J.,  1973:  Radar  observations  of  the  atmosphere.  Univ.  Of  Chicago  Press, 
Chicago,  HI. 

Bennetts,  D.  A.  and  B.  J.  Hoskins  1979:  Conditional  symmetric  instability  -  a  possible 
explanation  for  frontal  rainbands.  Quart.  J.  Roy.  Meteor.  Soc.,  105, 945-962. 

Boatman,  J.  F.  and  R.  F.  Reinking,  1984:  Synoptic  and  mesoscale  circulations  and 
precipitation  mechanisms  in  shallow  upslope  stonns  over  the  western  high 
plains.  Mon.  Wea.  Rev.,  112,  1725-1744. 

Btrfuie,  A.,  1979:  A  dual-Doppler  radar  study  of  mesoscale  winter  snow  bands.  Tech 
Repent  40, 196  pp.  Available  from  the  Laboratory  fen-  Atmospheric  Probing. 
Dept  .of  the  Geophysical  Sciences,  Univ.  of  Chicago,  Chicago,  Ill. 

Browning,  K.  A.  and  R.  Wexler,  1968:  The  detemnnation  of  kinematic  properties  of  a 
wind  field  using  Doppler  radar.  J.  Appl.  Meteor.,  7,  105-1 13. 


142 


Byrd,  G.  P.,  1989:  A  composite  analysis  of  winter  season  ovemining  precipitation 

bands  over  the  southern  plains  of  the  United  States.  J.  Atmos.  ScL,  46,  1 1 19- 
1132. 

Carbone,  R.  E.  and  A.  Bohne,  197S:  Cellular  snow  generation  -  A  Doppler  radar  study. 
J.  Attnos.  Sci.,  39,  1384-1394. 

Corby,  G.  A.  and  C.  E.  Wallington,  1956:  Airflow  over  mountains:  The  lee  wave 
amplitude.  Quart.  J.  Roy.  Meteor.  Soc.,  82,  266-274. 

Davies-Jones,  R.  P.,  1979:  Dual-Doppler  radar  coverage  area  as  a  function  of 

measurement  accuracy  and  spatial  resolution,  J.  Appl.  Meteor.,  18, 1229-1233. 

desJardins,  M.  L.,  R.  F.  Brill,  and  S.S.  Schotz,  1991:  GEMPAK5  User's  guide. 
NASA  Tech  Memo  4260. 

desJardins,  M.  L.,  and  R.  L.  Peterson,  1985:  GEMPAK:  A  meteorological  system  for 
research  and  education.  Preprints,  Int'l  Conf  on  Interactive  Information  and 
Processing  Systems  for  Meteorology,  Oceanography,  and  Hydrology,  Los 
Angeles,  American  Meteorological  Society,  313-319. 

Dirks,  R.  A.,  J.  P.  Kuetmer,  and  J.  A.  Moore,  1988:  Genesis  of  Atlantic  Lows 

Experiment  (GALE):  An  overview.  Bull.  Amer.  Meteor.  Soc.,  69,  148-172. 

Doviak,  R.  J.  and  D.  S.  Zmic,  1984:  Doppler  radar  and  weather  observations. 
Academic  Press  Inc.,  New  York,  458  pp. 

Duiui,  L.,  1987:  Cold-air  damming  by  the  Front  Range  of  the  Colorado  Rockies  and  its 
relationship  to  locally  heavy  snows.  Wea.  Forecasting,  2, 177-189. 

Durm,  L.  1988:  Vertical  motion  evaluation  of  a  Colorado  snowstcxm  from  a 
synoptician's  perspective.  Wea.  Forecasting,  3, 261-272. 

Durran,  D.  R.,  1986:  Mountain  Waves.  Chapter  20  in  Mesoscale  Meteorology  and 
Forecasting .  P.  S.  Ray,  Ed.  American  Meteorological  Society,  Boston,  MA 
472-492. 


143 


Emanuel,  K.  A.,  1983a;  On  assessing  local  conditional  symmetric  instability  from 
atmospheric  soundings.  Mon.  Wea.  Rev.,  Ill,  2016-2033. 

Emanuel,  K.  A.,  1983b:  The  Lagrangian  parcel  dynamics  of  moist  symmetric 
instability.  /.  Atmos.  ScL,  40,  2368-2376. 

Emanuel,  K.  A.,  1988:  Observational  evidence  of  slantwise  convective  adjustment 
Mon.  Wea.  Rev.,  116,  1805-1816. 

Hane,  C.  E.,  1986:  Extratropical  squall  lines  and  rainbands.  Chapter  16  of  Mesoscale 
Meteorology  and  Forecasting,  P.  S.  Ray,  Ed.  American  MeteOTological  Society, 
Boston,  MA  359-389. 

Heckman,  B.  E.  and  T.  Dulong  1989:  Doppler  radar  reflectivity  bands:  Comparison 
with  500  mb  patterns  and  implications  for  winter  nowcasting.  Preprints  3rd 
Cortf.  on  the  Aviation  Weather  System,  Anaheim,  CA.  233-231. 

Herzegh,  P.  H.  and  P.  V.  Hobbs,  1980:  The  mesoscale  and  microscale  structure  and 
organization  of  clouds  and  precipitation  in  midlatitude  cyclones.  Part  11: 
Warm-frontal  clouds.  J.  Atmos.  Sci.,  37, 597-61 1. 

Heymsfield,  G.  NL,  1979:  Doppler  radar  study  of  a  warm  frontal  region.  J.  Atmos. 
Sci.,  36. 2093-2107. 

Hobbs,  P.  V.,  T.  J.  Matejka,  P.  H.  Herzegh,  J.  D.  Locatelli,  and  R.  A.  Houze,  Jr., 
1980:  The  mesoscale  and  microscale  structure  and  organization  of  clouds  and 
precipitation  in  midiatitude  cyclones.  Part  I;  A  case  study  of  a  cold  front  J. 
Atmos.  Sci.,  37, 568-596. 

Hobbs,  P.  V.and  P.  O.  G.  Persson,  1982:  The  mesoscale  and  microscale  structure  and 
organizaticm  of  clouds  and  precipitation  in  midlatitude  cyclones.  Part  V:  The 
substructure  of  narrow  cold-frontal  rain  bands,  J.  Atmos.  Sci.,  39,  280-295. 

Hobbs  P.  V.  and  A.  L.  Rangno,  1985:  Ice  particle  concentrations  in  clouds.  J.  Atmos. 
Sci.,  42,  2523-2549. 


144 


Houze,  R.  A.,  S.A.  Rutledge,  T.  J.  Matejka,  and  P.  V.  Hobbs,  1981;  The  mesoscale 
and  microscale  structure  and  organization  of  clouds  and  precipitation  in 
midlatitude  cyclones.  Part  III;  Air  motions  and  precipitation  growth  in  a  warm- 
fiontal  rainband.  J.  Atmos.  Sci.,  38,  639-649. 

Kessinger,  C.  J.  and  W.  C.  Lee,  1991;  Evaluation  of  real-time  dual  Doppler  analysis 
for  use  during  field  operations.  Preprints  25th  Int'l.  Conf.  on  Radar 
Meteorology,  Paris,  France  756-759. 

Keyser,  D.  and  M.  A.  Shapiro,  1986;  A  review  of  the  structure  and  dynamics  of  upper- 
level  frontal  zones.  Mon.  Wea.  Rev.,  114,  452-499. 

Koch,  S£.,  M.  desJardin,  and  P.  J.  Kocin,  1983;  An  interactive  Barnes  objective  map 
analysis  scheme  for  use  with  satellite  and  conventional  data.  J.  Climate  Appl. 
Meteor.,  22,  1487-1503. 

Lhermitte,  R.  M  and  D.  Atlas,  1961;  Precipitation  motion  by  pulse  Doppler  radar. 

Proc.  Ninth  Weather  Radar  Conf.,  Boston,  MA.  Amer.  Meteor.  Soc.,  218-223. 

Lilly,  D.  K.,  1986;  Instabilities.  Chapter  1 1  of  Mesoscale  Meteorology  and  Forecasting, 
P.  S.  Ray,  Ed.  American  Meteorological  Society,  Boston,  MA  259-271. 

Lilly,  D.  K.,  1981;  Doppler  radar  observations  of  upslope  snowstorms.  Preprints,  2nd 
Corf,  on  Mountain  Meteorology,  AMS,  Steamboat  Springs,  CO.  Nov  1981, 
346-353. 

Lindzen,  R.  S.  and  K.  K.  Tung,  1976;  Banded  convective  activity  and  ducted  gravity 
waves.  Mon.  Wea.  Rev.,  104,  1^2-1617. 

Lcmg,  A.  B.,  B.  A.  Campistron,  and  A.  W.  Huggins,  1990:  Investigations  of  a  winter 
mountain  storm  in  Utah.  Part  I;  Synoptic  analyses,  mesoscale  kinematics  and 
water  release  rates.  J.  Atmos.  Sci.,  47,  1302-1322. 

Marwitz,  J.,  1980;  Winter  Storms  over  the  San  Juan  Mountains.  Part  I:  Dynamical 
processes.  J.  Appl  Meteor.,  19  913-926. 


145 


Marwitz,  J.  and  D.  Day,  1991:  The  effect  of  melting  in  a  Denver  snow  dump.  Preprints, 
First  Int'l  Winter  Storms  Symposium,  New  Orleans,  14  -  18  Jan  1991, 241- 
244. 

Matejka,  T.  J.,  R.  A.  Houze,  and  P.  V.  Hobbs,  1980:  Microphysics  and  dynamics  of 
clouds  associated  with  mesoscale  rainbands  in  extratropical  cyclones.  Quart.  J. 
Roy.  Meteor.  Soc.,  106,  29-56. 

Matejka,  T.  J.,  and  R.  C.  Srivastava,  1991:  An  improved  version  of  the  extended 

velocity-aamuth  display  analysis  of  single-Doppler  radar  data.  J.  Amos.  Ocean. 
Tech.,  8,  453-466. 

Mohr,  C.,  L.J.  Miller,  R.  Vaughan,  and  H.  Frank,  1986:  The  merger  of  mesoscale 

datasets  into  a  common  Cartesian  format  for  efficient  and  systematic  analysis.  J. 
Atmos.  Oceanic.  Technol.,  3,  143-161. 

Moore,  J.  T.,  and  P.  D.  Blakely,  1988:  The  role  of  ffontogenetical  fcH-cing  and 

conditional  symmetric  instability  in  the  midwest  snowstorm  of  30-31  January 
1982.  Mon.  Wea.  Rev.,  116,  2155-2171. 

Mulvey,  G.  J.,  1977:  Physical  mechanisms  of  extra  area  effects  fftHn  weather 

modification.  Phd  Dissertation,  Deptof  Atmospheric  Science,  Colraado  State 
University,  Fort  Collins,  CO,  138  pp. 

NCAR,  1991: 1991  Winter  icing  and  storms  project  (WISP)  data  catalog,  15  January  - 
5  April  1991,  Available  from  National  Center  for  Atmospheric  Research, 
Research  Applications  Program  Office,  Boulder,  CO. 

I^ffsons,  D.B  and  P.  V.  Hobbs,  1983:  The  mesoscale  and  microscale  structure  and 
Qi;ganization  of  clouds  and  precipitation  in  midlatitude  cyclones.  XI: 
Comparisons  between  observational  and  theoretical  aspects  of  rainbands.  J. 
Atmos.  Sci.,  40,  2377-2397. 


146 


Queney,  P.,  G.  Corby,  N.  Gerbier,  H.  Koschmieder,  and  J.  Zierep,  1960:  The  airflow 
over  mountains.  WMO  Tech.  Note  34, 135  pp. 

Raddatz,  R.  L.  and  M.  L.  Khandekar,  1979:  Upslope  enhanced  extreme  rainfall  events 
over  the  Canadian  western  plain:  A  mesoscale  numerical  simulation.  Mon.  Wea. 
Rev.,  107,  650-661. 

Ramamurthy,  M.  K.,  R.  M.  Rauber,  B.  P.  Collins,  M.  T.  Shields,  P.  C.  Kennedy,  and 
W.C.  Qark,  1991:  UNIWIPP:  A  University  of  Illinois  field  experiment  to 
investigate  the  structure  of  mesoscale  precipitation  in  winter  stcmns.  Bull.  Amer. 
Meteor.  Soc.,  72,  764-776. 

Raman,  S.  and  A.  J.  Riordan,  1988:  The  Genesis  of  Atlantic  Lows  Experiment:  The 

planetary-boundary-layer  subprogram  of  GALE.  Bull.  Amer.  Meteor.  Soc.,  69, 
161-172. 

Rasmussen,  R.  M.  and  M.  K.  Politovich,  1990:  WISP  Scientific  Overview.  NCAR 
publication  available  from  authors,  P.  O.  Box  3000,  Boulder,  Co  80307 

Rasmussen,  R.,  M  Politovich,  J.  Marwitz,  W.  Sand,  J.  McGinley,  J.  Snook,  R. 

Pielke,  S.  Rutledge,  D.  Wesley,  G.  Stossmeister,  B.  Bernstein,  K.  Elmore,  N. 
Powell,  E.  Westwater,  B.  Stankov,  D.  Burrows,  1992:  Winter  Icing  and 
Stcnms  Project.  Accepted  for  publication  in  Bull.  Amer.  Meteor.  Soc. 

Rauber,  R.  M.,  L.  O.  Grant,  D.  Feng,  and  J.  B  Snider,  1986:  The  characteristics  and 
distributions  of  cloud  water  over  the  mountains  of  northern  Colorado  during 
wintertime  storms.  Part  I:  Temporal  variations.  J.  Clim  Appl.  Meteor.,  25, 468- 
488. 

Reinking,  R.  F.  and  J.  F.  Boatman,  1986:  Upslope  precipitation  events.  Chapter  19  of 
Mesoscale  Meteorology  and  Forecasting,  P.  S.  Ray,  Ed.  American 
Meteorological  Society,  Boston,  MA  437-471. 


147 


Reuter,  G.  W.  and  M.  K.  Yau,  1990:  Observations  of  slantwise  convective  instability 
in  winter  cyclones.  Mon.  Wea.  Rev.,  118,  447-458. 

Reynolds,  D.  W.  and  A.  S.  Dennis,  1986:  A  review  of  the  Sierra  Cooperative  Pilot 
Project,  Bull.  Amer.  Met.  Soc.,  67,  513-523. 

Rioidan,  A.  J.,  1990:  Examination  of  the  mesoscale  features  of  the  GALE  coastal  front 
of  24-25  January  1986.  Mon.  Wea.  Rev.,  118,  258-282. 

Rogers,  D.  C.  and  G.  Vali,  1987:  Ice  crystal  production  by  mountain  surfaces.  J.  Clint. 
Appl.  Meteor.,  16,  1152-1168. 

Rutledge,  S.  A.  and  P.  V.  Hobbs,  1983:  The  mesoscale  and  microscale  strucmre  and 
otganization  of  clouds  and  precipitation  in  midlatitude  cyclones.  Part  Vni:  A 
model  for  the  "seeder-feeder"  process  in  warm-frontal  rainbands,  J.  Atmos. 
Sci.,  40,  1185-1206. 

Rutledge,  S.  A.  and  P.  V.  Hobbs,  1984:  The  mesoscale  and  microscale  structure  and 
organization  of  clouds  and  precipitation  in  midlatitude  cyclones.  Part  XII:  A 
diagnostic  modeling  study  of  precipitation  development  in  narrow  cold-frontal 
rainbands.  J.  Atmos.  Sci.,  41,  2949-2972. 

Sanders,  F.,  1986:  Frontogenesis  and  symmetric  stability  in  a  major  New  England 
snowstorm.  Mon.  Wea.  Rev.,  114,  1847-1862. 

Sanders,  F.  and  L.  F.  Bosart  1985a :  Mesoscale  structure  in  the  mega-lopolitan 
snowstorm  of  1 1-12  February  1983.  Part  I:  Frontogenetical  forcing  and 
symmetric  instability.  J.  Atmos.  Sci.,  42,  1050-1061. 

Sanders,  F.  and  L.  F.  Bosart  1985b :  Mesoscale  structure  in  the  mega-lopolitan 

snowstorm  of  1 1-12  February  1983.  Part  II:  Doppler  radar  study  of  the  New 
England  snowband  J.  Atmos.  Sci.,  42,  1398-1407. 

Sassen,  K.,  A.  W.  Huggens,  A.  B.  Long,  J.  B.  Snider,  and  R.  J.  Meitin,  1990: 

Investigations  of  a  winter  mountain  storm  in  Utah.  Part  II.  Mesoscale  structure, 


148 


supercooled  liquid  water  development  and  precipitation  processes.  J.  Atmos. 
Sci.,  47,  1323-1350. 

Schlatter,  T.  W.,  D.  V.  Baker,  and  J.  F.  Henz,  1983:  Profiling  Colorado's  Christmas 
Eve  blizzard.  Wecaherwise,  36, 60-66. 

Scorer,  R.,  1949:  Theory  of  waves  in  the  lee  of  mountains.  Quart.  J.  Roy.  Meteor. 
Soc.,  75,  41-56. 

Shapiro,  M.  A.,  1984:  Meteorological  tower  measurements  of  a  surface  cold  front 
Mon.  Wea.  Rev.,  112,  1634-1639. 

Shields,  M.T.,  R.  M.  Rauber,  and  M.  K.  Ramamurthy,  1991:  Dynamical  forcing  and 
mesoscale  organization  of  precipitation  bands  in  a  midwest  winter  cyclonic 
storm.  Mon.  Wea.  Rev.,  119,  936-964. 

Snook,  J.  S.,  1992:  Current  techniques  for  real-time  evaluation  of  conditional 

symmetric  instability.  Accepted  for  publication  in  Weather  and  Forecasting. 

Stivastava,  R.  C.,  T.  J.  Matejka,  T.  J.  Lorello,  1986:  Doppler  radar  study  of  the  trailing 
anvil  region  associated  with  a  squall  line.  /.  Atmos.  Sci,  43,  356-377. 

Stull,  R.  B.,  1988:  An  introduction  to  boundary  layer  meteorology.  Kluwer  Academic 
Publishers,  London. 

Szoke,  E.  J.,  1991 :  The  use  of  Doppler  radar  in  the  short-range  forecasting  of  snowfall 
in  northeastern  Colorado.  Preprints  25th  Int’l.  Conf.  on  Radar  Meteorology, 
Paris,  France  71-74. 

Todi,  J.  J.,  1991:  Surface  horizontal  pressure  gradients  in  upslope  winter 

stopaasJ*reprints  First  Int'l  Winter  Storms  Symporium,  New  Orleans,  14-18 
Jan.,  101-102. 

Toth,  J.  J.,  1987:  Interaction  of  shallow  cold  surges  with  topography  on  scales  of  100- 
1000  kilometers.  CIRA  Publication,  Colorado  State  University,  Fort  Collins, 
ISSN  No.  07373-5352-8,  135  pp. 


149 


Uccellini,  L.  W.  and  P.  J.  Kocin,  1987:  The  interaction  of  jetstreak  circulations  during 
heavy  snow  events  along  the  East  Coast  of  the  United  States.  Wea.  Forecasting, 
2,  289-308. 

Walsh,  P.  A.,  1977:  Cloud  measurements  in  wintertime  clouds.  M.  S.  thesis.  Dept  of 
Atmospheric  Science,  Univ.  of  Wyoming,  Laramie,  170  pp. 

Wesley,  D.  A.,  1991:  An  investigation  of  the  effects  of  topography  on  Colorado  Front 
Range  winter  storms.  Phd  Dissertation,  Dept.of  Atmospheric  Science,  Colorado 
State  University,  Fort  Collins,  CO,  201  pp. 

Wesley,  D.  A.  and  R.  A.  Pielke,  1990:  Observations  of  blocking-induced  convergence 
zones  and  effects  on  precipitation  in  complex  terrain.  Atmospheric  Research, 
25,235-276. 

Wesley,  D.  A.,  J.  F.  Weaver,  and  R.  A.  Pielke,  1990:  Heavy  snowfall  during  an 

extreme  arctic  outbreak  along  the  Colorado  Front  Range.  Natl.  Wea  Dig.,  15, 2- 
19. 

Whiteman,  C.  D.,  1973:  Some  climatological  characteristics  of  seedable  upslope  cloud 
systems  in  the  High  Plains.  NOAA  Tech.  Rep.  268  -APCL-27  (NTIS-COM-73- 
50924/2GI),  43  pp. 

Wolfsberg,  D.  G.,  K.  A.  Emanuel,  and  R.  E.  Passarelli,  1986:  Band  formation  in  a 
New  England  winter  storm.  Mon.  Wea.  Rev.,  114,  1552-1569. 


150 


APPENDIX  A 


Station  Identification  and  Location  Table 


151 


APPEMDIX  A 


Station  Idantlflcation  and  Location 

Balov  ara  tha  station  callsigns.  World  Meteorological  Organization 
station  nunbars  (if  assignad),  station  name,  state  Identification, 
country  identification,  latitude,  longitude  (-  indicates  west 
longitude)  and  alavatlon  in  meters.  Radar  sites,  wind  profilers,  and 
CLASS  stations  are  indicated  after  the  station  name.  All  stations, 
except  CHL  and  MHR  had  recorded  surface  weather  data. 
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T.aiMr 

CO 
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CO 

US 
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CO 

US 
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CO 

US 
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ARV 

99999 
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CO 

US 
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CO 

US 
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US 
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Briggsdala  PROFS 

CO 

US 

4065  -10434 

1480 
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Broomfield  (Jefferson  Co  AP) 

CO 

US 

3990  -10512 
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Aurora  (Buckley  AP) 

CO 
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3972  -10478 
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BOU 

99999 

Boulder  PROFS 

CO 
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Brighton  PROFS 
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US 
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US 
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Winter  Park 

CO 

US 
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Craig  (-Moffat  AP) 

CO 

US 

4052  -10755 
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MM 

US 

3645  -10315 
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US 
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WY 

US 

4115  -10482 
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CO 

US 

3975  -10487 
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CO 

US 

3965  -10691 

1985 
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99999 

Elbert  PROFS 

CO 

US 

3923  -10463 
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99999 
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CO 

US 
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CO 

US 
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US 

4058  -10508 
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US 
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US 
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US 
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US 
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US 
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Idaho  Springs  PROFS 
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US 
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3455 
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99999 

Klstell  Mobile  CLASS 

NE 

US 

4122  -10467 

1488 

xro 

99999 

Keansburg  PROFS 

CO 

US 

4007  -10451 

1519 

LAK 

99999 

Lakewood  PROFS 

CO 

US 

3970  -10516 

1825 

LAR 
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Laramie  (Gen  Braes  Fid) 
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US 

4132  -10568 

2217 

UF 

725620 

Worth  Platte  (Lee  Bird  Fid) 

NE 

US 

4113  -10068 

849 

LGM 

99999 

Longmont  PROFS 

CO 

US 

4017  -10516 

1536 

LHX 

724635 

La  Junta  (Muni  AP) 

CO 

US 

3805  -10353 
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LZC 
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CO 

US 

3930  -10370 
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LTN 

99999 

Littleton  PROFS 

CO 

US 

3957  -10496 
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LVB 

99999 

Loveland  PROFS 

CO 

US 

4041  -10504 

1513 

MHR 

99999 

HOAA  MHR  Radar 

CO 

US 

3987  -10478 
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mm 

99999 

Honn  PROFS 

CO 

US 

4080  -10476 

1638 

POl 

99999 

Tie  Siding  PAM 

WY 

US 

4106  -10554 
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P02 

99999 

Bums  PAM 

WY 

US 

4118  -10436 

1668 

P03 

99999 

Kimball  PAM/CLASS 

ME 

US 

4119  -10367 

1494 

P04 

99999 

Sidney  PAM 

NE 

US 

4110  -10298 

1307 

P05 

99999 

Rlst  Canyon  PAM 

CO 

US 

4061  -10530 

2210 

POt 

99999 

Owl  Canyon/Lew  Grant  Farm  PAM 

CO 

US 

4077  -10509 

1657 

P07 

99999 

Carr  PAM 

CO 

US 

4090  -10488 

1743 
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P08 

99999 

H«r«£ord/Poly  B«rger  PAM 

CO 

US 

4097  -10429 

1615 

P09 

99999 

Battl*  Canyon  PAM 

CO 

US 

4095  -10368 

1526 

PIO 

99999 

Stonaham  PAM 

CO 

US 

4070  -10369 

1361 

Pll 

99999 

Starling  PAM 

CO 

US 

4061  -10321 

1202 

P12 

99999 

Barttaoud  PAM/CIASS 

CO 

US 

4032  -10514 

1577 

P13 

99999 

Mastars  PAM 

CO 

US 

4029  -10424 

1369 

P14 

99999 

Higgins  PAM/CIASS 

CO 

US 

4029  -10405 

1384 

P15 

99999 

Midway  PAM 

CO 

US 

4022  -10340 

1356 

P16 

99999 

Akron  PAM/CLASS 

CO 

US 

4017  -10321 

1420 

P17 

99999 

Nadarland  PAM 

CO 

US 

3995  -10551 

2589 

P18 

99999 

Battasso/Bouldar  Filtar  PAM 

CO 

US 

4001  -10534 

1951 

P19 

99999 

North  of  Bannatt  PAM 

CO 

US 

3991  -10440 

1564 

P20 

99999 

Hoyt  PAM 

CO 

US 

4002  -10413 

1460 

P21 

99999 

Noodrow  PAM 

CO 

US 

3999  -10357 

1359 

P22 

99999 

Staplaton  PAM 

CO 

US 

3977  -10487 

1606 

P23 

99999 

Front  Ranga  Airport  PAM 

CO 

US 

3979  -10456 

1672 

P24 

99999 

South  of  Bannatt  PAM 

CO 

US 

3956  -10445 

1795 

P25 

99999 

Last  Chanca/Lusto  SP  PAM 

CO 

US 

3961  -10358 

1570 

P26 

99999 

Anton  PAM 

CO 

US 

3972  -10312 

1457 

P27 

99999 

Sadalla  Damon  Butta  PAM 

CO 

US 

3937  -10497 

1881 

P28 

99999 

Bllzabatb  PAM 

CO¬ 

US 

3943  -10462 

1981 

P29 

99999 

Bijou  Saum  Ranch  PAM 

CO 

US 

3941  -10424 

1852 

P30 

99999 

Agata  PAM 

CO 

US 

3946  -10394 

1664 

P31 

99999 

MonusMnt  PAM 

CO 

US 

3909  -10486 

2146 

P32 

99999 

Blbart  PAM/CLASS 

CO 

US 

3923  -10463 

2135 

P33 

99999 

Mathason  PAM 

CO 

US 

3905  -10398 

1874 

P34 

99999 

Lljaon  PAM 

CO 

US 

3918  -10370 

1694 

P35 

99999 

Flaglar  PAM/CLASS 

CO 

US 

3936  -10304 

1480 

P36 

99999 

Bllcott  PAM 

CO 

US 

3887  -10441 

1876 

P37 

99999 

Truckton/Todar  PAM 

CO 

US 

3877  -10417 

1803 

P38 

99999 

Punkln  Cantar  PAM 

CO 

US 

3884  -10388 

1692 

P39 

99999 

Hlld  Horsa  PAM 

CO 

US 

3891  -10298 

1408 

PLV 

99999 

Plattavllla  Moblla  CLASS 

CO 

US 

4022  -10480 

1515 

PTL 

99999 

Plattavllla  PROFS/Hlnd  Profllar 

CO 

US 

4026  -10487 

1449 

PUB 

724640 

Puablo  (Mast  AP) 

CO 

US 

3828  -10451 

1439 

RAT 

99999 

Haw  Rayisar  Moblla  CLASS 

CO 

US 

4061  -10384 

1430 

ROL 

99999 

Rolllnsvllla  PROFS 

CO 

US 

3991  -10549 

2748 

RNL 

725745 

Rawlins  (Muni  AP) 

HY 

US 

4180  -10720 

2055 

SMT 

725610 

Sldnay  (Muni  AP) 

NB 

US 

4110  -10298 

1312 

TAD 

724645 

Trinidad 

CO 

US 

3725  -10433 

1750 

WU> 

99999 

Hard  PROFS 

CO 

US 

4004  -10554 

3004 

III 

99999 

99999  99999 

99999 
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APPENDIX  B 


Summary  of  Data  Display  Development 

The  GEneral  Meteorological  data  assimilation,  analysis  and  display  software 
PAcKage  (GEMPAK)  developed  by  the  National  Aeronautics  and  Space  Administration 
(NASA)  was  used  extensively  in  the  analysis  of  meteorological  data  particularly, 
conventional  surface  and  upper  air  meteorological  data  (desJardins  and  Peterson,  1985; 
Koch  et  al.  1983).  The  breadth  and  utility  of  GEMPAK  programs  saved  cmisiderable 
time  in  analysis  since  virtually  any  data  type  or  parameter  can  be  analyzed  and  displayed 
with  (xily  minor  modifications  to  the  GEMPAK  command  files  (desJardins  er  oZ.  1991). 
Most  of  the  figures  in  this  thesis  showing  conventional  data  were  prepared  using 
GEMPAK.  This  required  writing  FORTRAN  programs  to  convert  a  wide  variety  of 
Hata  foiciats  to  GEMPAK  format  as  well  as  initializing  and  organizing  GEMPAK  data 
files.  The  primary  data  source  was  NCAR.  Their  WISP  91  Data  Catalog  (NCAR; 
1991)  gives  a  ccmqtlete  description  of  data  collected  during  the  experiment  Conversion 
programs  were  written  for  the  following  data  sets: 

PortaUe  Automated  Mescmet  (PAM)  surface  observations 
Prototype  Regional  Operational  Forecast  System  (PROFS)  surface  observations 
Gross-chain  Loran  Atmospheric  Sounding  System  (CLASS)  rawinsonde 
observations 

National  Weather  Service  rawinsonde  observations 
30  second  U.  S.  Geological  Survey  topographic  data 
Wind  profiler  observations 
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Data  was  also  analyzed  and  displayed  using  the  PLOT  GKS  graphics  package 
develq)ed  at  Oregon  State  University  and  later  modified  and  expanded  at  Colorado  State 
University.  Time  -  height  series.  Scorer  parameter,  stability  analyses,  and  some  radar 
data  displays  woe  plotted  with  this  program.  FORTRAN  programs  were  written  to  put 
the  the  fcdlowing  into  network  common  data  format  -  netcdf  files  used  by  the  PLOT 
GKS  program: 

PAM  pressure  data 
Profiler  data 

30  second  U.  S.  Geological  Survey  topographic  data 
Extoided  Vertical  Azimuth  Display  (EVAD)  data 
Time-height  cross  secticxis  of  EVAD  derived  divergence,  vertical  velocity, 
Drqjpler  radial  velocity,  and  reflectivity  fields 
Dual  Dcqipler  analyses  of  vertical  velocity,  horizontal  velocity,  divergence,  and 
reflectivity  fields 

In  additicm  to  PLOT  GKS,  other  programs  were  used  to  edit  and  display  radar 
dafp  (4)tained  fiom  the  CSU  CHILL  radar  and  the  Mile  High  Radar  (MHR).  Raw  data 
were  checked  for  fcdded  Doppler  velocities  and  unfolded  if  necessary  interactively  using 
the  NCAR  developed  Research  Data  Support  System  (RDSS)  editing  software  on 
UNDC*  based  SUN*  wcnkstations.  Vertical  cross  sections  were  also  prepared  and 
examined  using  tiiis  software.  After  ctxnpletion  of  editing,  data  were  converted  from 
radar  qtherical  coordinates  to  Cartesian  cocxdinates  with  the  Soned  Position  Radar 
INTeipc^laticm  (SI^INT)  software  package.  Finally,  the  Custom  Editing  and  Display 
of  Reduced  Information  in  the  Cartesian  Space  (CEDRIC)  Doppler  radar  analysis 
jsogram  was  used  for  single-  and  dual-Doppler  data  analysis  and  creation  of  netcdf  data 
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files.  See  Mohr  et  ai.  (1986)  for  more  information  on  SPRINT  and  CEDRIC.  The 
netcdf  files  generated  by  SPRINT  and  CEDRIC  were  called  by  PLXDT  GKS  for  data 
display.  Control  of  the  displays  was  accomplished  through  manipulation  of  a  command 
file  called  on  execution  of  PLOT  GKS. 

RDSS  was  used  in  conjunction  with  the  X' Windows  XWD,  San  Diego 
Supercomputer  Center  IMCONV,  and  Wavefront  IMF_DSPL  programs,  on  a  Silicon 
Graphics  Iris®  Personal  Work  Station  to  generate  loops  of  MHR  reflectivity  and  velocity 
fields  at  the  Colorado  State  University  Visualization  Laboratory.  The  loops  were  then 
wiitren  to  a  video  disk  for  subsequent  recording  on  videotape.  Also,  an  overlay  of 
terrain  contours  obtained  fiom  the  30  second  U.  S.  Geological  Survey  terrain  height  data 
was  included  on  the  loc^.  This  provided  insight  into  location  and  movement  of  echoes 
with  respect  to  major  terrain  features. 
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APPENDIX  C 

Colorado  State  University  CHILL  and  Mile  High  Radar  Specifications 


CSU'CHBLL  Radar  Characteristics  in  WISPQl 


System: 

Doppler  CapaUlity  yes 

Systdn  Gain  (dB)  (Assumed  to  be  antenna  gain  minus  wareguide 
loss.  Also  CHILL  virtually  never  operates  w/o  the  radome] 
ZDR  Mode  39 

Bypass  Mode  41 

Polarisation  [can  vary  from  pulse  to  pulse] 

Xmit  H  Rec  H 
Xmit  VRecV 


Antenna: 

Shape 

Fhed 

Diameter  (m) 

Beamwidth  (H/V  deg) 
Antenna  Gain  (dB) 

First  Sidelobe  (1  way  db) 
Scan  Rate  (deg/s) 

Anmuth 

Elevation 

Noise  Temperature  (k) 

(indttding  waveguide] 

Tkansmitter: 

Wavdength  (cm) 
Fkequenqr  (MBs) 

Peak  Power  (kW) 

Pulse  Width  (micro-scc) 
PRF(Hs) 


parabolic 

center 

8.5 

.96 

43.3  dB 

-25  [radome  on] 

6-24 

0-15 

»100 


11 

2725 

600 

1.0 

1050  {WISP91  average] 
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CSU'CHILL  Radar  Charaeteriities  -  continued 


Receiver: 

Noise  temp  (K)  »  420 

Ndse  Figure  (dB)  4.0 

Trensfer  function 

Doppler  linear 

Intensity  linear 

Log  Channel  [Did  not  exist  in  WISP91] 

Linear  channel 

Bandwidth  (MBs)  5 

Dynamic  Range  (dB)  50 

Noise  Power  (dBm)  -111 .  1 

Min  detectable  Signal  (dBm)  »  —125 

Min  Det  Refi.  Factor 

[rain  O  50  km,  dBZ  -15  [based  on  S/N  of  -15  dB] 

Data  Acquisition: 

No.  of  range  gates 
Range  gate  spacing  (m) 

No.  of  Samples 

Recorded  (R)  and  Displayed  (D)  Data  Fields: 


Log  Power  hor. 

R,D 

Normalised  Coherent  Power 

R,D 

Spectrum  Width 

R,D 

Mean  radial  velocity 

R,D 

500-1400 

150 

55-128  [range  for  WISPOl] 
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Capabilities 


MHR 

System 

Doppler  Capability 

yes 

System  Gain  (dB) 

Radome  on 

43 

Radome  off 

— 

Polarization 

Xmt  Rec  H 

Antenna 

Shape 

parabolic 

Feed 

center 

Diameter  (m) 

8.2 

Beamwidth  (H/V  deg) 

0.9 

Antenna  gain 

44 

First  Sidelobe  (1  way  dB) 

Radome  on 

-25 

Radome  off 

Scan  rate  (deg/s) 

Azimuth 

0-30 

Elevation 

0-15 

Noise  temperature  (K) 

«100 

(includes  waveguide) 

Transmitter 

Wavelength  (cm) 

10 

FVequency  (MHz) 

Peak  Power  (kW) 

700 

Pulse  Width  (us) 

1.6 

PRF  (Hz) 

250-1250 

Receiver 

Noise  temp.  (K) 

400 

Noise  figure  (dB) 

3.8 

Transfer  function 

Doppler 

IF  linear 

Intensity 

IF  linear 

Log  Channel 

Bandwidth  (MHz) 

— 

Dynamic  Range  (dB) 

— 
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MHR 

Noise  power  (dBm) 

Linear  Channel 

— 

Bandwidth  (MHz) 

0.6 

Dynamic  Range  (dB) 

100 

Noise  power  (dBm) 

-107 

Min.  Det.  Signal  (dBm) 

-122 

Min.  Det.  Refl.  Factor 
(Rain  0  50  km,  dBz) 

-18 

Data  Acquisition 

No.  of  range  gates 
(intervals  of  256  gates  ) 

256-1024 

Rajage  gate  spacing  (m) 

225-900 

No.  of  samples 

45, 14 

Recorded  (R)  and  Displayed  (D)  Data  Fields 

Log  Power,  hor. 

Coherent  power,  hor. 

R,  D 

Spectrum  vddth 

R,D 

Velocity 

Log  Power,  circ.  main  power 

Log  Power,  circ.  cross  polar. 

Coherent  Power,  circ.  main 

CDR 

Correlation 

R,  D 
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